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Abstract 


A  detailed  review  of  the  seven  lightning  return 
stroke  models  is  presented.  These  models  were  proposed  by 
Bruce  and  Golde  (1964),  Dennis  and  Pierce  (1964),  Little 
(1978),  Y.T.  Lin  et  al.  (1980),  R.L.  Gardner  (1980), 

Strawe  et  al.  (1980),  and  Master  et  al .  (1981).  For  the 
better  understanding  of  the  subject,  a  review  of  lightning 
return  stroke  process  is  also  included. The  equations 
to  compute  the  electric  and  magnetic  fields  in  space  due 
to  a  vertical  channel  are  also  derived.  The  last  part 
of  the  thesis  is  related  to  the  comparison  of  these  models. 
The  predicted  fields  due  to  Master  et  al.  were  plotted 
and  compared  with  the  empirical  results.  Als’o,y‘ Master 
et  al.  model  (1981)  was  modified  for  the  channel  current 
pulse  propagation  velocity  and  the  channel  current  wave 
form.  The  plots  due  to  these  modifications  are  included 
and  compared  with  the  experimental  data. 


A  REVIEW  AND  COMPARISON  OF  LIGHTNING 
RETURN  STROKE  MODELS  USING 
EXPERIMENTAL  DATA 

I.  Introduction 


«r» 


Over  the  last  three  decades  many  lightning  return 
stroke  models  have  been  proposed  to  predict  the  current  pulse 
of  the  return  stroke  as  it  propagates  in  the  channel.  These 
models  were  mainly  proposed  by  trying  to  solve  the  inverse 
problem  of  suggesting  the  current  in  the  channel  which  will 
give  the  experimental  measurement  of  the  electromagnetic  (EM) 
fields  on  the  ground.  Due  to  the  lack  of  airborne  data,  the 
validity  of  these  models  as  a  function  of  height  have  not 
been  verified.  In  this  thesis  we  provide  .^i  discussion  of  all 
the  significant  return  stroke  models  and  their  credibility  as  a 
function  of  height  by  comparing  their  EM  fields  with  recent 
acquired  airborne  lightning  data.'(Rustan  et  al.,  1982). 

Prior  to  1970,  the  main  concern  about  the  lightning 
phenomena  was  due  to  its  effect  on  transmission  lines.  How¬ 
ever,  from  1970  to  1981  seven  USAF  aircraft  were  lost  due  to 
confirmed  lightning  related  incidents  and  more  than  150 
serious  lightning  related  mishaps  were  recorded  (Corbin: 

IEEE  Trans  1982).  These  numerous  incidents  coupled  with 
recent  changes  of  the  aircraft  structure  to  composite  ma¬ 
terials  alarmed  the  electromagnetic  interference  community 
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of  the  ever  present  lightning  hazards.  Additional  testing  of 
the  aircraft  to  the  airborne  lightning  environments  was  re¬ 
quired.  The  new  fiber-reinforced  plastics  and  other-  non- 
metallic  materials,  which  have  replaced  the  conventional  alum¬ 
inum  structure,  are  more  susceptible/vulnerable  to  lightning 
damage  because  these  are  usually  incapable  of  withstanding 
high  voltage  stresses  created  by  lightning  and  thus  causing 
dielectric  breakdown  (J.A.  Plumer  &  J.D.  Robb,  1982),  Although 
the  use  of  composite  materials  in  aircraft  technology  has 
resulted  in  lower  cost,  lighter  weight  and  lower  radar  return 
signatures,  it  has  also  caused  undesirable  external  EM  energy 
coupling  into  the  aircraft.  Particularly,  the  external  EM 
radiation  at  high  frequency  is  a  serious  hazard  for  modern 
electronic  devices,  which  are  more  responsive  at  higher  fre¬ 
quencies  and  can  be  destroyed  at  lower  voltages  than  their 
predecessors  (Corbin, 1979) .  Electromagnetic  transients  pro¬ 
duced  by  lightning  currents  have  frequency  spectrum  from 
near  DC  to  GHz  range  (Uman  et  al:1980).  These  transients 
when  coupled  inside  the  aircraft  produce  resonance  effects, 
particularly  in  the  range  of  1  MHZ  to  10  MHZ.  This  in  turn 
critically  affects  the  performance  of  wideband  digital  sys¬ 
tems  which  are  not  shielded  properly.  In  summary,  return 
stroke  EM  radiation  may  jeopardize  the  aircraft  safety  by 
causing  data  and  software  error  in  the  onboard  computers 
(Baum : 1980) . 

To  determine  the  EM  coupling  to  the  aircraft  for 
direct  and  nearby  lightning,  one  must  have  an  accurate  model 


of  the  lightning  process.  By  doing  so  one  can  determine  the 
level  of  interference  and  what  phase  of  the  lightning  discharge 
is  likely  to  produce  induced  electromagnetic  pulses.  This  will 
help  in  designing  the  aircraft  structure  and  shielding  for  the 
modern  digital  electronics  systems  inside  the  aircraft. 

In  order  to  determine  the  needed  protection  and  to  learn 
more  about  the  physical  properties  of  the  lighting  discharge, 
many  researchers  have  developed  lightning  return  stroke  models 
which  describe  different  aspects  of  this  natural  phenomena. 
Within  the  scope  of  this  thesis,  the  author  has  chosen  to  re¬ 
view  the  most  significant  work  in  this  area.  The  return  stroke 
models  to  be  included  in  our  discussions  are  the  ones  proposed 
by 


(a) 

C.E.R.  Bruce  and  R.H.  Golde 

(1941) 

(b) 

A.S.  Dennis  and  E.T.  Pierce 

(1964) 

(c) 

P.F.  Little  (1978) 

(d) 

Y.T.  Lin,  M.A.  Uman  and  R.B. 

S^randler  (1980) 

(e) 

R.L.  Gardner  (1980) 

(f) 

P.F.  Strawe  et  al  (1980) 

Lightning  return  stroke  models  derived  by  the  above 
mentioned  authors  predict  EM  field  at  ground  level.  However, 
M.J.  Master  et  al  (1981)  proposed  a  model  to  predict  EM  fields 
at  altitudes.  This  model,  in  fact,  is  a  slightly  modified 
version  of  Lin's  model,  which  he  presented  in  his  Ph.D. 


thesis. 


Problem 


The  problem  of  this  thesis  is  three-fold.  First  of 
all,  the  most  significant  lightning  return  stroke  models  found 
in  the  literature  are  reviewed  and  critically  analyzed. 
Secondly,  the  results  from  the  Master's  model  are  analyzed 
and  compared  with  the  flight  data  obtained  near  active 
thunderstorms  at  the  same  altitudes  and  distances.  This  will 
enable  us  to  determine  whether  or  not  the  proposed  model  re¬ 
produces  the  surface  fields  measured  on  an  aircraft  during 
airborne  tests.  Finally,  we  will  modify  the  Master's  model 
in  the  form  of  variations  in  certain  parameters  to  establish 
relationship  between  lightning  currents  and  the  resultant  EM 
fields.  Since  most  of  the  previous  models  use  data  measured 
at  ground  level  for  their  analysis,  we  will  determine  which 
one  of  these  models  best  reproduces  the  electromagnetic  fields 
as  a  function  of  height. 

Scope 

The  scope  of  this  thesis  is  to  collect  and  fully 
analyze  the  information  contained  in  different  lightning 
return  stroke  models  in  concise  and  compact  form.  This  will 
enable  the  researcher  Interested  in  lightning  phenomena  to 
have  quick  access  to  some  of  the  important  models  presented 
over  the  last  few  decades.  The  analysis  and  comparison  part 
of  this  thesis  will  enable  the  reader  to  select  a  model  for 
further  studies  and  modifications.  In  selecting  the  return 
stroke  lightning  models  which  appear  most  significant,  we  have 
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taken  equal  contributions  of  the  linear  and  the  nonlinear 
approaches  into  model  development. 

Lin  et  al  (1980)  model,  with  modification  by  Master 
et  al  (1981),  is  presently  the  only  model  which  derives  EM 
fields  as  a  function  of  height  and  distances.  This  model  is 
studied  in  detail  in  this  thesis.  The  comparison  of  their 
results  with  the  data  collected  by  the  Air  Force  Flight  Dyna¬ 
mic  Laboratory  will  show  how  well  this  model  can  predict  the 
fields  produced  by  a  lightning  return  stroke. 

The  Masters  et  al,  computer  model,  used  to  predict  the 
fields  produced  by  a  lightning  current  as  a  function  of  height 
and  distances  is  based  on  the  statistical  averages  of  the 
model  parameters  described  by  Lin  et  al  (1980).  Additional 
variations  in  the  channel  current  parameters  are  introduced 
to  calculate  the  fields.  By  comparing  the  results  of  these 
variations  with  the  experimental  data  we  can  determine  the 
refinements  required  in  this  model.  Due  to^the  time  constraint; 
in  this  study  it  was  not  possible  for  the  author  to  explore 
all  the  possible  variations  and  adjustments  in  this  model. 

Presentation 

In  Chapter  II  the  lightning  flash  mechanism  is  descri¬ 
bed.  A  study  of  return  stroke  model  essentially  requires  the 
understanding  of  this  chapter.  Relevant  terminologies  used 
in  lightning  literature  are  also  defined  here.  Chapters  III 
and  IV  describe  the  linear  and  the  nonlinear  models,  respec¬ 
tively.  Equations  for  the  calculation  of  EM  field  at  a  point 


in  space  are  derived  in  Chapter  V.  These  derivations  will 
provide  an  analytic  insight  to  the  Master  et  al.  model. 

Master  et  al.  model  (1982)  and  modification  to  it  are  discussed 
in  Chapter  VI  and  Chapter  VII  discusses  the  validity  of  experi¬ 
mental  data.  Comparison  of  models  is  presented  in  Chapter  VIII. 
Finally,  Chapter  IX  summarizes  the  conclusions  and  recommenda¬ 
tions  based  on  this  research. 


II.  Lightning  Mechanism 


Cloud-to-Ground  Lightning 

Lightning  is  a  transient,  high  current  electric  dis¬ 
charge,  which  has  a  path  length  extended  over  several  kilo¬ 
meters  (Uman,  1969).  It  is  commonly  initiated  by  thunder¬ 
clouds,  which  range  from  one  or  two  km  diameter  to  giant  elec¬ 
trical  storms  as  shown  in  Fig  2.1  (Vonnegut,  1965).  Most  of 
the  lightning  flashes  that  occur  around  the  world  are  intra¬ 
cloud  discharges  (70  to  80^  of  all  lightning  in  Florida),  but 
the  type  of  lightning  flashes  which  are  of  most  concern  to 
us  are  cloud-to-ground  and  these  are  the  only  ones  studied. 

Generally,  cloud-to-ground  lightning  is  referred  to 
as  streaked  or  forked  lightning.  An  entire  discharge  to 
ground  lasts  typically  .2  to  .25  seconds  -and  is  called  a  flash 
(Bruce  Gold,  1941;  Uman,  1969).  The  flash  is  composed  of  in¬ 
dividual  discharges  called  strokes,  each  stroke  lasts  for  a 
few  milliseconds  with  a  time  between  strokes  of  tens  of 
milliseconds.  There  are  typically  three  or  four  strokes  per 
flash,  the  minimum  being  one,  the  maximum  measured  26.  Fig  2.2 
shows  a  distribution  of  the  number  of  strokes  in  a  flash  in 
different  parts  of  the  world.  Thomson  (1980)  reports  that 
there  is  no  significant  correlation  between  the  average  number 
of  strokes  per  flash  and  the  average  interstroke  time  interval. 
There  are  varying  parameters  which  affect  the  lightning 
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Fig  2.1  Comparison  of  Various  Sizes  of  Convective  Clouds 
that  Produce  Lightning  Discharge  (adapted  from 
Vannegut,  1965). 
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Fig  2.2  Number  of  Successive  Strokes  in  Lightning  Flashes 
(Adapted  from  Bruce  &  Golde,  1941). 


;  characteristics.  Frontal  storms,  for  example,  produce  a  higher 

flash  rate  and  more  strokes  in  a  flash  than  local  convective 
storms  (Uman,  1982).  Similarly,  the  relative  height  of  the 
terrain  affects  the  channel  length. 

In  1930  Schonland  et  al. ,  using  photographic  techniques, 
determined  that  lightning  flashes  were  initiated  by  the  stepped 
leaders  in  the  channel.  To  explain  this  phenomena  a  simple 
model  for  the  charge  structure  of  a  thunderstorm  was  made  from 
ground-based  measurements  of  the  cloud  electric  fields  and 
lightning  field  changes.  Fig  2.3  shows  such  a  model,  in  which 
the  major  region  of  positive  charge  build-up  occurs  at  the 
top  of  the  cloud  (P) .  Below  this  positive  charge  distribution 
there  is  a  strong  negative  charge  distribution  (N).  The 
O’  bottom  of  the  cloud  has  a  weak  positive  charge  distribution 

(p) .  Although  the  actual  process  by  which  stepped  leader 
propagate  is  unknown,  the  local  electrical  breakdown  between 
the  N  and  p  region  of  the  thundercloud  is  cpnsidered  to  cause 
this  phenomena.  Streak-camera  photographs  as  shown  in  Fig  2.4 
indicate  that  the  stepped  leader  has  a  strongly  luminous 
step.  There  is  a  silent  period  between  steps  and  each  new 
step  has  a  faintly  luminous  channel  in  the  clouds.  Thus,  the 
stepped  leader  is  a  cloud-to-ground  predischarge  preceding 
the  first  return  stroke. 

The  average  length  of  a  stepped  leader  is  50  meters 
with  pause  time  of  50  microseconds  (ysec)  between  steps. 

^  The  average  velocity  of  steps  towards  the  ground  is  nearly 

•  e: 

1.5x10  meters/sec.  However,  Baum  et  al  (1980)  suggested  that 
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Fig  2 


- - - - - 

(•)  (b) 


4  (a)  Streak  Camera  Photograph  of  Stepped  Leader 

Process,  Scale  of  Drawing  is  Distorted 
for  Illustrative  Purpose. 

(b)  Ordinary  Camera  Picture  of  a  Lightning 
Flash.  (Adapted  from  Uman,  1969). 
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stepped  leaders  may  form  in  shorter  than  10  m  steps  and 
faster  than  the  velocity  calculated  by  earlier  authors.  When 
stepped  leaders  propagate  down  towards  earth  they  lower  a 

Q 

negative  potential  of  the  order  of  -10  volts  and  on  the  aver¬ 
age  they  distribute  negative  charge  of  5  coulombs  along  their 
path. 

When  the  stepped  leader  nears  the  ground,  the  electric 
field  there  is  sufficient  to  cause  streamers  to  be  initiated 
and  propagate  upward  to  meet  the  leader  tip.  At  the  junction 
of  the  downward-moving  leader  and  the  upward  moving  streamer 
the  "return  stroke"  wavefront  forms  and  thereafter  propagates 
upwards  carrying  ground  potential  into  the  cloud.  The  return 
stroke  is  responsible  for  ground  and  channel  currents  of 
typically  20  kA  rising  to  peak  value  in  about  one  microsecond. 
However,  risetirae  to  peak  current  of  hundreds  of  nanoseconds 
has  also  been  reported  in  the  literature  (Berger  et  al.  1975). 
The  above  measurements  are  at  the  base  of  the  channel,  measure- 
ments  of  the  channel  current  above  the  ground  are  lacking 
(Price  &  Pierce  1977). 

The  upward  velocity  of  the  return  stroke  is  of  the 
8 

order  of  1x10  meter/sec  and  the  trip  between  ground  and 
cloud  base  is  completed  in  approximately  70  usee.  Bruce  and 
Golde  (1941)  suggested  that  the  velocity  of  the  return  stroke 
decreases  as  it  proceeds  upwards  and  it  can  be  represented 
as 

V  =  Vq  e  (2.1) 
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8  4 

where  ,8x10  meters/sec  and  y  ~  3x10  .  (If  t  =70  psec 

then  according  to  this  relation  the  velocity  of  the  return 
stroke  at  the  bottom  of  the  cloud  is  9.796x10°  m/sec.) 

Lin  (1978)  also  reports  that  the  return  stroke  velocity  near 

Q 

the  ground  is  1x10  meters/sec  and  at  the  top  of  the  channel 
is  of  the  order  of  4x10^  meters/sec.  Boyle  and  Orville 
(1976)  suggested  that  this  velocity  is  of  the  order  of  .05  - 
.5  of  the  speed  of  light.  This  variation  in  velocity  with 
height  indicates  that  there  is  possible  dispersion  and  dis¬ 
sipation  of  the  return  stroke  pulse  along  the  channel. 

As  the  charge  in  the  channel  is  depleted,  the  stroke 
current  ceases  to  flow  ending  the  flash.  However,  if  the 
channel  is  replenished  with  charge  (from  the  N  region  in  clouds ) 
within  about  100  microseconds,  a  continuous  or  "dart  leader” 
may  start  moving  down  the  previous  return  stroke  channel. 

The  downward  velocity  of  the  dart  leader  is  almost  uniform 

g 

(of  the  order  of  2x10  m/sec)  and  greater  than  that  of  the 
stepped  leader.  The  return  stroke  caused  by  the  dart  leader  is 
termed  as  "a  subsequent  return  stroke." 

Sometimes  the  lower  part  of  the  stroke  channel  dis¬ 
appears  due  to  long  interstroke  periods  or  high  wind  veloci¬ 
ties  shear  and  disintegrate  the  channel.  In  such  cases 
the  dart  leader  on  its  downward  trip  changes  its  role  to  a 
stepped  leader  form.  This  kind  of  leader  is  named  as  "dart 
stepped  leader”  and  it  initiates  the  first  return  stroke. 


Subsequent  return  stroke  field  changes,  though  of 
smaller  magnitude,  are  similar  to  but  faster  than  those  of 
the  first  return  stroke.  Also  the  subsequent  return  stroke 
has  faster  current  risetime  and  similar  maximum  rate  of 
change  of  current.  Consecutive  return  strokes  in  the  same 
channel  are  usually  separated  by  40  to  80  microseconds; 
however,  this  interval  may  exceed  100  msec  provided  a  continuous 
current  keeps  on  flowing  in  the  channel  (Uman  &  Krider,  1981). 

We  have  briefly  described  the  characteristics  of  a 
return  stroke  which  lowers  negative  potential  and  negative 
charge  from  the  cloud  towards  ground.  The  return  stroke 
models  discussed  in  this  thesis  are  of  this  form.  Occasionally, 
stepped  leaders  lower  positive  charge.  Return  strokes  caused 
by  this  type  of  leader  have  slow  rate  of  rise  of  current  and 
large  charge  transfer.  Positive  discharge  flashes  are  rarely 
composed  of  more  than  one  stroke.  Another  kind  of  stepped 
leader  is  that  which  moves  upward  towards,. the  cloud  and  may 
carry  either  positive  or  negative  charge.  A  detailed  descrip¬ 
tion  of  this  form  of  stroke  is  given  by  Uman  (1969). 

Return  Stroke  Current  Waveforms 

In  the  modeling  of  return  strokes,  the  current  waveform 
has  a  key  role.  As  said  earlier,  the  direct  current  measure¬ 
ments  have  always  been  done  at  the  ground  level.  The  most 
common  device  to  measure  the  amount  of  current  and  its  wave¬ 
form  is  the  magnetic  link.  This  device  is  preferred  because 
it  can  measure  peak  currents  regardless  of  their  duration  and 


relative  variations.  However,  the  fast  risetime  to  peak 
currents  or  fast  variations  in  the  current  waveform  are  limited 
by  the  system  response. 

In  the  early  works  on  lightning,  the  peak  value  of  the 
first  return  stroke  current  was  measured  as  15-20  kA.  (McCann, 
1944;  Lewis  8s  Foster,  1945).  The  accuracy  of  their  measure¬ 
ments  was  questioned  by  Anderson  and  Hegenguth  (1959)  and 
Szpor  (1969).  Szpor  (1969)  suggested  that  these  values  were 
low  by  a  factor  of  two  and  thus  peak  current  value  is  approxi¬ 
mately  30  kA  at  50%  distribution  level.  Similar  values  for 
peak  currents  were  calculated  by  Berger  etal^,(1975)  and 
Garbagnati  et  al .  ,(1974) .  The  average  peak  value  in  subse¬ 
quent  return  strokes  is  almost  half  of  that  in  the  first  re¬ 
turn  stroke. 

Bruce  and  Golde  (1941)  suggested  that  the  return  stroke 
waveform  could  be  represented  as 

I(t)  =  (e"“^  -  e"^^)  (2.2) 

where  a  =  4.4x10^  ,  6  =  4.6x10^,  =  30  kA  and  I(t)  is  the 

current  at  the  base  of  the  return  stroke  at  time  t.  These 
values  give  average  time  to  crest  and  to  half  value  of  about 
6  psec  and  23.5  psec , .respectively .  Current  waveforms  observed 
by  different  authors  differ  by  various  parameters  and  such 
waveforms  are  shown  in  Fig  2.5.  Parameters,  which  are  compared 
to  these  waveforms,  are  risetime  to  peak  current,  time  to 
half  value  of  peak  current,  and  charge  transferred.  The 
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( b )  SUBSEQUENT  STROKE  CURRENT  WAVEFORM 


Fig  2.5  Typical  First  and  Subsequent  Return  Stroke 

Current  Waveforms:  Averaged  Current  by  Berger 
et  al.  (1975).  Typical  Currents  by  Hagenguth 
et  al.  (1952)  and  Garbagnatic  et  al.  (1974). 
(Adapted  from  Lin,  1978). 


current  waveform  of  a  lightning  return  stroke  at  the  base 
of  the  channel  has  been  suggested  by  Uman  (1969)  as 


I  =  ^o 


_  .-St  ^ 


^1  ® 


(2.3) 


This  waveform  differs  from  that  of  the  Bruce  8i  Golde 
(1941)  by  the  quantity  which  is  considered  as  inter¬ 

mediate  current.  For  the  first  return  stroke  values  of  1^, 

4—1  5-1 

CL  and  6  are  30  kA,  2.0x10  sec  and  2.0x10  sec  ,  respectively 

In  a  subsequent  return  stroke  these  values  are  10  kA,  1.4x 
4  -1  4  -1 

10  sec  and  6.0x10  sec  .  The  intermediate  current  has  been 

3  -1 

specified  as  =  2.5  kA  and  y  =  1.0x10  sec 

From  Fig  2.4  and  datalisted  above,  it  is  apparent  that 
the  first  return  stroke  has  longer  risetime,  higher  peak  cur¬ 
rent,  longer  half-value  time  for  current  waveform  and  larger 
amount  of  charge  transferred  than  the  subsequent  return  strokes 
There  are  other  parameters  in  the  first  return  stroke  which 
differ  from  the  subsequent  return  stroke,  such  as  channel 
tortuousity  and  velocity  of  propagation.  In  this  thesis  the 
channel  will  be  considered  as  straight  and  vertical  whereas 
velocity  parameter  will  be  discussed  in  one  of  the  next 
chapters . 
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III.  Linear  Models 


This  chapter  analyzes  the  linear  return-stroke  models 
proposed  by  Bruce  &  Golde  (1941),  Dennis  &  Pierce  (1964)  and 
Y.T.  Lin  et  al . ,  (1980).  These  models  provide  a  relationship 
between  return-stroke  EM  fields  and  the  corresponding  channel 
current.  These  models  attempt  to  solve  the  inverse  problem  of 
estimating  the  current  in  the  channel  based  on  the  measurements 
of  the  EM  fields  a  certain  distance  away  from  the  channel. 

If  the  estimated  current  in  the  channel  satisfies  other  mea¬ 
surements  of  the  EM  fields  as  a  function  of  distance  we  claim 
the  model  is  adequate.  Each  one  of  these  three  authors  pur¬ 
sued  this  same  objective  and  kept  making  improvements  over 
previous  models. 

Bruce  and  Golde  Model  (1941) 

In  this  model  the  return  stroke  current  at  any  given 
time  is  assumed  to  be  uniform  with  height  below  the  return 
stroke  wave  front  and  zero  above  it.  Thus,  the  current  at 
the  base  of  the  channel  and  below  the  wave  front  are  identical, 
that  is 


and 


i(z,t)  =>  i(o,t)  z  £  1  (3.1) 

i(z,t)  =  o  z  ^  1  (3.2) 
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where 


i(z,t)  =  IqJcxp  (-at)  -  exp  (-Bt)] 


(3.3) 


is  the  current  wave  form  at  the  ground,  and  is  the  height 
of  the  channel.  Also,  Bruce  and  Golde  (BG)  assumed  from  the 
photographic  data  of  Schonland  (1956)  that  for  the  first  re¬ 
turn  stroke  the  propagation  velocity  decreases  with  time  and 
is  given  by 


Vt  =  Vq  exp  (-yt)  (3.4) 

(Values  of  y  and  v^  are  given  in  Chapter  II.). 

However,  for  the  subsequent  return  stroke  this  velocity 
was  assumed  to  be  constant.  An  illustration  of  Bruce-Golde  model 
is  given  in  Fig  3.1. 

To  calculate  fields  from  the  above  current  wave  form 
Bruce-Golde  used  the  "charge  moment"  equations  (Lejay,  1926; 
Ratcliffe  et  al.,  1932).  According  to  these  equations  the 
electric  and  magnetic  field  strength  at  a  distance  R  from  a 
dipole  of  moment  M  are  given  by 


E^CR.t) 


fM 

.  1 

.  1 

d^Ml 

cR^ 

dt 

“5” 

c^R 

dt^J 

B  (R,t) 


^o  r  1  dM  ^  1  d^Ml 


(3.5) 


(3.6) 


where  "c"  is  the  velocity  of  light,  is  the  permittivity 
of  the  medium  and  y  is  the  permeability  of  the  medium. 


Fig  3.1  Current  Distribution  in  Bruce  &  Golde  Model  (1941) 
(Adapted  from  Lin,  1978) 


In  equation  (3.5)  the  first  term  is  the  electrostatic  field,  the 
second  term  is  the  induced  field,  and  the  last  term  is  the 
radiated  field.  The  time  rate  of  change  of  electric  moment 
is  given  by 


dt 


2i(t) 


dt 


(3.7) 


The  factor  2  allows  for  image  charges.  Substituting  equations 
(3.3)  and  (3.4)  in  equation  (3.7),  yields 


^  =  2  J  exp  (-at)  -  exp  (-et)jj  1-exp  (-Yt)j  (3.8) 

Fields  derived  from  equations  (3.5)  and  (3.6)  in  con¬ 
junction  with  equation  (3.8)  agreed  relatively  well  with  the 
measured  values  (Lin,  1978). 

However,  there  are  two  obvious  deficiencies  in  this 
model.  First  of  all,  since  it  is  assumed  that  the  current 
distribution  along  the  channel  is  uniform,  the  charge  and  the 
magnitude  of  the  current  at  the  tip  of  the  channel  are  required 
to  be  instantaneously  transferred  from  the  base  of  the  channel, 
that  is,  infinite  velocity.  This  phenomena  contradicts  physi¬ 
cal  laws  because  the  travel  velocity  cannot  exceed  the  speed 
of  light.  Secondly,  there  is  always  a  current  discontinuity 
at  the  wave  front.  This  causes  an  infinite  rate  of  change  of 
current  for  zero  time.  The  second  problem  can  be  solved  by 
considering  the  finite  rate  of  change  of  the  current  and  then 
taking  the  limit  which  converges;  however,  the  first  problem 
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cannot  be  physically  justified.  In  order  to  correct  the  infi 
nite  speed  problem,  Dennis  Si  Pierce  (1964)  modified  the  Bruce 
Golde  model,  which  will  be  discussed  in  the.  following  section 

Dennis  and  Pierce  Model  (1964) 

From  equation  (3.4)  the  channel  length  z  at  time  't' 
is  given  by 

z  =  /  dt  (3.9) 

•b 

To  remove  the  discrepancy  in  the  BG  model,  Dennis  & 
Pierce  modified  equation  (3.3)  so  that 

i(z,t)  =  i(t-z/u)  (3.10) 

where  oi*  is  the  assumed  constant  velocity  of  current  waveform 
behind  the  return-stroke  wavefront,  and  th^  quantity  z/u  pro¬ 
vides  the  delay  time  in  getting  the  information.  This  model 

becomes  the  BG  model  if  is  set  equal  to  infinity.  Dennis 
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Si  Pierce  considered  values  of  u  of  8x10  ,  3x10  and  infinite 
meter /sec  in  their  computations.  However,  the  value  of  u  in 
their  final  model  was  equal  to  the  speed  of  light. 

Dennis  &  Pierce  used  charge  deposit  equations,  which 
will  be  defined  next,  to  calculate  the  electric  moment  M  and 
its  derivatives.  To  derive  charge  equation,  consider  that  at 
any  time  't'  the  charge  at  the  tip  of  the  channel  is  "frozen" 
or  deposited  there  to  neutralize  the  opposite  charge  stored 


earlier  at  that  point  by  the  stepped  leader.  Let  the  channel 
at  any  time,  t,  has  a  flowing  charge  density  p  per  unit  length 
and  deposited  charge  density  p'  per  unit  length;  then  the 
total  charge  per  unit  length  will  be  (p+p').  Assuming  that 
p  is  a  function  of  the  channel  height,  h,  and  time,  t,  and 
that  the  current  waveshape  proposed  by  BG  in  equation  (3.3) 
is  adequate,  then  the  charge  deposit  can  be  computed  by  the 
equation 


If  p^  is  the  charge  density  at  the  channel  tip  of 
height  z,  then  the  charge  deposited  at  the  tip  of  the  channel 
in  a  unit  time  is  whereas  the  charge  density  deposited 

in  a  unit  length  is 


where 


0^.  =  — 
t  y 


exp  [ -a(t-z/u)]  -  exp  [-B(t-z/u)j 


(3.13) 


From  equation  (3.9)  in  general 


Substituting  these  values  in  equation  (3.7)  gives 
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dM  =  2  i^  ^  dh 
t 

But  —  is  the  total  charge  (p+p*)  in  the  channel  of  unit  length 
t 

as  defined  earliers.  Hence 


dm 


2  (p+p')  h 


Electric  moment  M  due  to  the  total  charge  in  channel 
length  of  height  z  is  obtained  by  integrating  the  above  equa¬ 
tion.  Thus 

Ulr  z 

M  =  2  j  (p  +  p' )  h  dh  (3.14) 

o 

Since  p'  is  not  a  function  of  time  and  p,p'  are  zero 
at  z  =  0,  by  differentiating  the  above  equa.tion  using  Leibnitz's 
Rule,  we- obtained . 


dM  _ 
dt 


hdh  +  2 


/^(P+P')  ~  dh  +  2  (p^+p'^) 


i 


.dz 

^dt 


(3.15) 


In  addition,  as  the  channel  height,  h,  is  not  a  func¬ 
tion  of  time,  and  the  second  term  in  (3,15)  is  zero,  the  result 
is 
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(3.16) 


^  _ 

dt 


L 


dh 


+ 


2(p  +pV)z 

t  ^ 


dz 

dt 


Derivatives  of  equations  (3.9)  and  (3.11)  once  sub¬ 
stituted  in  equation ^3 . 1^  give 


2I„u 


2I.U 


^  exp(-at)  [exp(^)J  -1 - ^  exp(-6t) 


[exp(i|)  -l] 


(3.17) 


Further  differentiation  of  equation  (3.17)  yields 


^  =  2IqU  exp(-at)  [l-exp(^']+  2I^v^  exp(-Yt)  exp 

[-«(t-^)]  -2IqU  exp(-Bt)  [l-exp~)]  -2IqV^^  exp 

(-Yt)  exp  [-3(t-^)]  (3.18) 

Radiated  field  calculated  from  equation  (3.18)  at  a 
distance  of  100  km  is  shown  in  Fig  (3.2),  and  is  compared  with 

the  BG  model  in  which  u  =  “.  The  values  of  a,  8,  y  and  used 

-4  -1  5-14  -1 

for  computation  are  4.5x10  sec  ,  4.5x10  sec  ,  3x10  sec 

and  30  kA,  respectively. 

Fig  3.2  shows  that  the  peak  radiation  field  intensity 
decreases  with  decrease  in  velocity,  u,  and  also  the  peak  value 
is  delayed  in  time.  However,  the  radiation  fields  become  almost 
identical  after  100  psec . 


Since  the  maximum  length  of  a  return  stroke  channel  is 


(•I 

(Cl  U*1iI0’iii/mc 


■ 

1 

I 

■ 


<p 


't 

t 

k  • 


I  ■ 

[4 

r 

I 

k; 

b 

L  ■ 


•  ^ 

—  to  increase  the  length  either  v  has  to  be  increased  keeping 

Y 

Y  constant  or  y  is  to  be  decreased  keeping  constant.  Dennis 
&  Pierce  chose  the  latter  option  (third  possibility  of  varying 
Vq  and  Y  simultaneously  not  considered).  Fig  3,3  shows  that 

no  appreciable  effect  in  the  radiated  field  occurs  due  to  vari¬ 
ations  in  the  channel  height.  The  value  of  y  is  restricted  by 
equation  (3.4). 

The  Dennis  8c  Pierce  (DP)  model  has  also  zero  current 
above  the  wavefront  and  thus  a  discontinuity  which  was  considered 
a  shortcoming  in  the  BG  model.  Therefore  the  basic  physical 
limitations  of  the  BG  model  are  not  solved  by  the  DP  model. 

In  order  to  remove  this  problem  Uman  and  McLain  (1969)  pro¬ 
posed  a  return  stroke  model,  which  is  known  as  "Transmission 
Line  Model"  in  the  literature.  Since  this  transmission  line 
model  is  the  basis  for  Lin's  model  (1978),  it  will  be  briefly 
discussed  here. 

Transmission  Line  Model 

In  this  moiel  the  channel  current  takes  the  form 
i(z,t)  =  Ig^exp  [-a(t-|)]  -  exp  [-6(t-^) 

where  equation  (3.19)  is  of  the  same  form  as  equation  (3.10). 

The  difference  between  the  two  current  waveform  lies  in  the 
velocity  of  propagation  of  the  wavefront.  In  equation  (3.19) 
the  wavefront  travels  along  the  channel  with  a  velocity,  v, 
as  it  would  along  a  suitable  transmission  line,  Uman  and 
McLain  (1969)  used  arbitrary  current  waveshapes  in  their  model 
but  in  this  thesis  it  will  be  the  same  as  suggested  by  BG, 


)! 


(3.19) 
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Although  the  transmission  line  model  is  physically  reasonable 
and  provides  easy  analytical  computation  of  EM  fields,  the 
current  in  this  model  does  not  lower  the  charge  stored  in 
the  corona  sheath.  Instead,  the  model  represents  impulse¬ 
like  current  due  to  the  atmospheric  breakdown  between  the  end 
of  the  leader  and  start  of  the  return  stroke.  Fig  3.1  illus¬ 
trates  the  current  in  this  model,  whereas  Fig  3.4  provides 
the  geometry  used  to  derive  the  field  equations. 

Radiation  fields  at  a  distance  D  (D>50  kM)  from  the 
bottom  of  a  straight  vertical  channel  are  proportional  to  the 
product  of  the  return-stroke  velocity  and  its  current  wave¬ 
shape  (Uman  &  McLain  1969;  and  Uman  and  McLain  1971)  for  any 
time  't'  such  that  t<H/v,  where  H  is  the  height  of  channel. 
This  yields 


perfect  conducting  plane  earth  cannot  be  computed  at  any 
distance  by  the  use  of  these  equations.  For  return-stroke 
current  i(z,t)  the  fields  at  any  distance  can  be  calculated 
from  the  equations  derived  by  Oman  &  McLain  (1969  &  1971). 


E^CD 


^  }  2  3Sin\  ^  -R/c)dT  dz 


+  (  2-3Sin_9  ^  ^  t-R/c)  dz 

I  cR^ 


'J  fif  It 


(3.22) 


B^(D.t) 


H 


^  It  i 


(3.23) 


since  by  virtue  of  equation  (3.4) 


=/  ''t 


Therefore,  variable  velocity  of  propagation  can  be 
incorporated  in  these  equations.  Detailed  derivations  of 
equations  (3.22)  and  (3.23)  are  given  by  Uman  et  al.  (1975). 


31 


Lin's  Model  (198^ 


To  remove  the  shortcomings  of  the  Bruce  Golde,  Dennis 
and  Pierce  and  Transmission  Line  models,  a  new  model  was  pro¬ 
posed  by  Lin  et  al.  (1980).  The  channel  current  in  the  new 
model  is  composed  of  three  separate  components:  (1)  a  short- 
duration  upward  propagating  break  down  current  pulse  of  con¬ 
stant  magnitude  and  waveshape.  This  pulse  is  responsible  for 
the  return  stroke  peak  current  and  the  risetime.  This  pulse 
is  assumed  to  propagate  up  at  a  constant  velocity;  (2)  a  uni¬ 
form  current  which  may  already  be  flowing  (leader  current) 
or  it  may  start  to  flow  soon  after  the  return  stroke  begins; 

(3)  a  corona  current  caused  by  the  downward  movement  of  charge 
initially  stored  in  the  corona  sheath  around  the  leader  channel. 

The  current  contribution  due  to  the  three  components  is 
illustrated  in  Fig  3.5.  The  geometry  of  the  channel  is  the 
same  as  for  the  transmission  line  model  shown  in  Fig  3.4. 

The  channel  current  at  any  instant  time  't'  at  altitude 
z  and  due  to  source  at  z'  of  length  Az '  is  represented  by 


i(z,z' ,t)  =  i^+I 


exp(-a(t-t ' ) 


-exp(-6 (t-t ' )  11  Az' 


(3.24) 


for  t^t '  and  z'^z  where  i^  is  the  current  in  the  breakdown  pulse. 
The  time  t'  is  defined  by  the  relation 


f  =  +  t  + 

u  on  c 


(3.25) 
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Current  Distribution  in  Lin  et  al.,  Model  (1980) 
(Adapted  from  Lin  et  al.,  1980) 


where  is  the  time  in  which  the  pulse  current  i|^  attains 

its  maximum  value.  The  constant  in  Equation  (3.24)  has 

typical  values  of  10-50  A/meter.  Values  of  a  and  6  are 

5  —1  6  —1 

arbitrarily  set  equal  to  10  sec  and  3x10  sec  respectively, 

and  the  scale  factor  X  is  assumed  to  have  value  of  1200. 

In  this  mode  the  value  of  I^,  the  uniform  current,  is 

dE 

determined  from  the  slope  ^  of  the  electric  field  near  the 
lightning  during  the  linear  ramp  region  from  the  following 
equation : 


u' 


=  2it€, 


2  2^/2 
(H-^+D^)  dE 

3t 


(3.26) 


The  above  equation  should  be  used  with  some  reservation 
to  calculate  I^.  Uniform  current  field  and  electrostatic  field 
are  both  distance  dependent.  At  a  distant  station  the  electric 
field  at  100  ysec  due  to  becomes  very  small  and  thus 
the  ramp  of  the  field  is  difficult  to  detect.  A  calculation 
of  from  such  a  field  becomes  highly  erratic  and  overesti¬ 
mated  in  magnitude.  Also,  if  the  field  is  measured  much  closer 
to  the  lightning  channel  (distance  <  2  kM)  then  the  ramp  is 
more  effected  by  the  corona  current.  Measurement  in  such  a 
condition  will  also  be  inaccurate  and  overestimated. 

The  corona  current  waveshape  is  assumed  to  be  identical 
at  all  heights  but  its  magnitude  decreases  exponentially  with 
height.  Corona  current  in  the  channel  flows  towards  the  ground 
at  the  speed  of  light.  A  fast  turn  on  time  corona  current  is 
considered  the  best  for  this  model.  A  good  agreement  between 


predicted  and  measured  fields  is  obtained  when  the  corona 


current  reaches  its  peak  value  in  1  usee  after  being  turned 
on  by  the  breakdown  pulse  and  decays  exponentially  with  a 
10  ysec  time  constant.  It  may  be  noted  that  the  corona 
current  waveform  at  any  altitude  along  the  channel  is  the 
B.G.  current  waveform  which  decreases  exponentially  in  am¬ 
plitude  with  height.  Analytically,  this  waveform  is  repre¬ 
sented  by  the  last  part  of  equation  (3.24)  with  the  values  of 
a  and  6  mentioned  on  page  34. 

Once  the  current  in  the  vertical  lightning  channel  is 
established  then  the  associated  fields  can  be  computed  from 
Equations  (3.22)  and  (3.23).  However,  the  corona  current  at 
altitude  z,  due  to  the  charge  sources  above  it,  can  be  computed 
by  integrating  Equation  (3.24)  with  respect  to  z'  from  z  to  H. 

A  simplified  procedure  in  this  respect  is  proposed  by  Lin  et 
al.  (1980). 

This  model  provides  an  insight  to  the  channel  losses 
which  are  not  considered  in  lossless  transmission  line  model. 
The  EM  fields  predicted  in  this  model  are  much  closer  to  the 
experimental  data  than  those  predicted  using  the  previous 
models. 
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IV,  Nonlinear  Models 


The  linear  models  which  were  presented  in  the  last 
chapter  did  not  account  for  certain  important  aspects  of  the 
lightning  phenomena.  For  example,  most  of  the  lightning  chan¬ 
nels  are  tortuous  and  tortuosity  in  a  channel  is  a  random  fac¬ 
tor.  There  are  not  straightforward  techniques  that  can  be  used 
to  approximate  the  tortuosity  in  the  channel.  However,  the 
effect  of  tortuosity  on  the  radiated  EM  fields  can  be  estimated. 
Also,  in  the  previous  models  earth  was  considered  as  perfect 
conductor  whereas  it  has  finite  amount  of  conductivity.  Radia¬ 
tion  in  the  earth  ionosphere  environment  is  yet  another  factor 
which  has  not  been  considered  in  linear  models.  Robert  L. 
Gardner  (1980)  considers  all  these  parameters  in  a  nonlinear 
lightning  channel  model  consisting  of  arbitrarily  connected  and 
oriented  line  segments. 

To  account  for  time  varying  parameters  in  a  return 
stroke  channel,  P.F,  Little  (1978)  describes  a  nonuniform  tran- 
mission  line  model  using  a  lumped-parameter  network.  In  his 
paper  he  describes  methods  to  calculate  such  parameters.  These 
calculations  will  be  presented  in  brief  form  during  the  des¬ 
cription  of  this  model.  Variation  in  the  current  pulse  during 
return  stroke  propagation  is  also  reviewed  in  his  paper. 

The  third  nonlinear  transmission  line  model  reviewed 
in  this  thesis  was  presented  by  Strawe  et  al.  (1980).  The 


authors  of  this  model  predict  the  return  stroke  current 
rise  rate  and  velocity  of  propagation  on  the  basis  of  non¬ 
linear  breakdown  physics  of  the  channel  arc.  The  current 
waveform,  as  in  other  models,  is  not  based  upon  experimental 


data.  In  fact,  the  current  waveform  in  this  model  is  the 
solution  of  energy  balance  equations,  which  are  the  result 
of  physical  phenomena  occurring  in  an  arc  channel.  In  this 
model,  it  is  claimed  that  current  rise  rate  decays  rapidly 
with  propagation  distance  from  the  point  of  return  stroke 
initiation. 

For  the  ease  of  understanding,  the  first  model  re¬ 
viewed  in  this  chapter  is  by  P.F.  Little  (1978),  the  next 
model  is  of  R.L.  Gardner  (1980)  and  finally,  the  model 
presented  by  Strawe  et  al.  will  be  reviewed. 

Little's  Model  (1978) 

In  this  model,  the  lightning  chanpel  is  formed  of 
resistive  and  inductive  elements  in  a  series  RLC  circuit 
as  shown  in  Fig  4.1.  The  switch  SW  shown  in  this  figure 
connects  the  channel  to  ground  where  the  final  stepped 
leader  meets  the  short  return  leader  propagating  upward 
from  the  ground  level.  Thus,  a  charged  transmission  line 
is  terminated  in  a  resistance,  and  the  channel  charge  is 
allowed  to  flow  towards  the  ground.  It  is  assumed  that 
the  inductance  in  this  connecting  path  is  negligible  and 
hence,  a  fast-rising  current  pulse  is  produced.  However, 


the  magnitude  of  current  peak  reduces  as  the  channel  height 
increases  above  the  ground. 


Since  the  transmission  line  presented  in  Pig  4.1 
consists  of  nonuniform  lumped  RLC  elements,  the  shape  of 
the  current  pulse  will  change  with  height. 

To  find  the  voltage  and  current  at  the  ith  element 
of  the  RLC  network  shown  in  Fig  4.1,  consider  the  trans¬ 
mission  line  equations; 


and 


G 


(4.1) 


(4.2) 


Differentiating  equation  (4.1)  and  substituting 

(4.2)  yields 


a^v. 

— 

3Vi 

3^V^ 

3z^ 

—  KC 

3t 

3t^ 

a^ii 

3Ii 

TO 

—  KC 

3t 

+  llC 

3t^ 

(4.3) 


(4.4) 


The  solutions  to  equations  (4.3)  and  (4.4)  deter¬ 
mine  the  value  of  voltage  and  current  in  the  transmission 
line  at  height  z  and  time  t. 

P.F.  Little  uses  uniform  transmission  line  seg¬ 
ments.  In  that  case,  if  V  is  the  input  voltage  at  the  base 
of  the  channel,  which  is  terminated  in  its  characteristic 
impedance,  the  current  at  the  ground  rises  to  its  maximum 


value  depending  on  R.  Since  the  transmission  line  Is  uni¬ 
form,  the  subscript  In  equation  (4.3)  can  be  dropped 
and  it  takes  the  form 


ifv  _  pf,  av  , 


(4.5) 


The  RLC  parameters  are  derived  from  physical  argu¬ 
ments  rather  than  from  analytic  treatment.  To  determine 
the  capacitance  element,  consider  the  vertical  channel  as 
a  charged  conductor  of  small  radius  and  having  the  same  po¬ 
tential  all  along  the  channel.  Depending  upon  the  size  of 
the  clouds  and  Its  height  above  the  ground,  the  channel  can 
be  modeled  as  shown  In  Fig  4,2,  The  electrostatic  field 
distribution  Is  determined  by  numerical  method  and  the 
results  are  shown  in  Fig  4,3. 

By  determining  the  electrostatic  E-field  around  a 
cylindrical  channel  at  a  distance  r,  the  charge  on  the 
channel  per  unit  length  A  becomes 


X  =  2'nre^E 


(4.6) 


and  the  capacitance  per  meter,  C  ,  is  calculated  from  the 
relation 


C*  =  ^ 
^  V 


(4.7) 


The  tip  nearest  the  ground  can  be  considered 
spherical  shape  containing  the  charge 
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Fig  4.1  Lumped  Parameter  Transmission  Line  as  a  ,  Bepresentative 
of  a  Lightning  Return  Stroke  in  Little's  Model  (1978). 


(Q)  (b) 


Fig  4.2  Lightning  Channel  Models  for  Small  and  Large  Clouds 
(a)  Spherical  Cloud  Model  (b)  Plane  Cloud  Model 
(Adapted  from  Little,  19'^8), 


Radius  km 

(a) 


Fig  4. 


(b) 


3  Electrostatic  Field  Lines  Around  the  Lightning 
Channel  of  Fig  4.2  (a)  Spherical  Cloud 

(b)  Plane  Cloud  (Little,  1978) 


It  is  to  be  noted  that  the  charge  distribution  should 
affect  the  length  of  the  gap  to  be  covered  by  the  last  step¬ 
ped  leader  due  to  the  height  of  gap  above  the  ground  and  the 
dimensions  and  height  of  the  cloud.  Typical  charge  distribu¬ 
tions  in  the  channel  are  shown  in  Fig  4.4a  and  Fig  4.5b, 

Once  C  is  found  then  for  any  channel  segment  the 
associated  can  be  calculated  by  multiplying  C  with  the 
length  of  the  segment.  However,  the  value  of  changes  with 
time.  This  is  because  the  value  is  associated  with  the 
beginning  of  the  return  stroke  as  the  field  distribution 
changes,  the  value  of  also  changes  in  each  segment.  This 
argument  puts  this  model  in  the  category  of  nonlinear  models. 
However,  P.F,  Little  has  preferred  to  use  the  constant  values 
of  during  the  channel  discharge. 

The  displacement  current  between  the  ground  and  the 
line  segments  of  the  model  provides  a  closed  loop  path  for 
the  segment  current.  Assuming  that  the  field  distribution 
is  stationary  then  the  displacement  current  is  given  by 

Id  =  S  H 

Now  the  inductance  associated  with  this  current 

n 

flow  in  the  nth  segment  is  determined  by  the  area  under  the 


Fig  4.4a  Distribution  of  Linear  Charge  Density  X  Along 
3  km  Lightning  Channel  for  Two  Clouds  Model 
(Little,  1978). 


Fig  4.4b  Distribution  of  Charge  Density  X  Along  3  km  and 
6  km  Lightning  Channels  with  Gap  Length  of  100  m 
at  an  Altitude  of  100  m  (Little  1978). 


field  lines.  The  inductance  for  closed  path  is  given  by 

1  7-  - 

L  =  ^qB.ds. 

where  is  considered  as  constant.  Assuming  that  the 
channel  radius  is  very  small  as  compared  to  the  segment 
length,  the  area  S  can  be  approximated  by  an  arc  of  radius 
equal  to  the  segment  length  as  shown  in  Fig  4.5.  The  shaded 
and  dotted  zones  in  Fig  4.5  determine  S  for  the  nth  element 
at  the  beginning  of  the  return  stroke.  Inductance  due  to 
(n+l)th  element  is  calculated  in  the  same  manner  and  associ¬ 
ated  area  is  the  dotted  region  of  Fig  4.5.  The  difference 
between  the  two  areas,  shown  as  hatched  portion  in  Fig  4.5a 
then  determines  the  flux  due  to  nth  element.  Circular  or 
rectangular  approximation  within  the  same  limits  as  shown 
in  Fig  4.5a,  b  give  essentially  the  same  value  of  L^,  The 
inductance  of  line  segment  is  then  determined  by  con¬ 
sidering  the  channel  as  coaxial  line,  whose  outer  radius 
is  equal  to  the  length  of  the  segment.  P.F.  Little  (1978) 
approximates  the  radius  of  the  channel  equal  to  2  cm  and 
assuming  that  field  lines  are  circular,  he  suggests  that 
inductance  per  unit  length  is  roughly  2  yH/m. 

The  resistance  element  in  the  channel  is  difficult 
to  calculate  because  it  is  a  complex  function  of  time  and 
height.  R*,  the  resistance  per  meter,  is  estimated  in  such 

a  manner  that  the  current  remains  unidirectional  but  its 

* 

oscillations  are  not  completely  dampened.  A  value  of  R 
of  1  ohm  per  meter  is  assumed  in  this  model.  In  the  switch 
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Gap  length  & 
height  above 
.ground  25  m. 

Gap  length  & 
height  above 


Time  (|is) 


ground  100  m. 


in  a  3  km  Lightning  Channel 
Height  Above  Ground  (m):  U 
(c)  1950,  (d)  1450,  (e)  900. 
(h)  75,  (Little,  1978) 
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Fig  4.6b  Current  Waveforms  in  a  3  km  Lightning  Channel 
for  a  Plane  Cloud  Model.  Height  Above  Ground 
(m):  (a)  2950,  (b)  1950,  (c)  900,  (d)  450, 

(e)  225,  (f)  75.  Gap  length  100  m  and  200-300  m 
Above  Ground  (Little,  1978). 


Fig  4.6c  Current  Waveforms  in  3  km  Lightning  Channel  for 
a  Plane  Cloud  Model.  Gap  Length  of  100  m  and 
600-700  ra  Above  Ground.  Height  Above  Ground  (m): 
(a)  2950,  (b)  1950,  (c)  825,  (d)  675,  (e)  450, 

(f)  150.  (Little,  1978). 
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loop  the  resistance  is  approximated  as  600  ohm.  This 
value  represents  the  energy  loss  in  the  breakdovm  process. 

Finally,  the  terminating  resistance  is  given  by 

Rg  =  8.9  (4.10) 

max 

where  p  is  the  specific  resistivity  of  the  ground. 

Equation  (4.10)  gives  a  value  of  R„  =  100  for  = 

£i  luaX 

-  100  kA.  The  cloud  capacitance  is  assumed  to  be  1  pF. 

Even  though  the  paraimeter  values  predicted  by  the 
above  model  are  not  correct,  they  are  fairly  close  to  the 
values  observed  on  the  ground.  Considering  this  aspect  of 
the  model,  the  current  pulse  shape  of  the  channel  current 
can  be  approximated  at  all  altitudes.  The  computed  current 
waveforms  at  different  heights  are  shown  in  Fig  4.6a,  b  &  c. 


Gardner's  Model  (1980) 

This  model  considers  the  lightning  channel  as  a  num¬ 
ber  of  arbitrarily  oriented  line  segments’  connected  together 
Contrary  to  other  models  which  consider  the  earth  as  a  per¬ 
fect  conductor,  this  model  accounts  for  the  finite  conduc¬ 
tivity  of  the  earth.  Also,  the  EM  fields  radiation  in  the 
earth  ionosphere  waveguide  are  included  in  this  model.  The 
current  in  the  channel  is  treated  as  a  pulse  when  a  pre¬ 
charged  transmission  line  is  subjected  to  a  discharge.  The 
current  pulse  pareuneters  are  the  seune  as  described  in 
equation  (2.3) .  The  geometry  of  a  point  charge  on  the 


channel  is  shown  in  Fig  4.7. 

In  order  to  calculate  the  E  and  5  fields  due  to  a 

line  segment,  the  three  orthogonal  components  L^,  L^,  and 

are  calculated  where  L  ,  L„  and  represents  the  length  of 

X  y  z 

the  segment  component  in  x,  y,  and  z  direction.  The  current 
pulse  is  assumed  to  travel  along  these  orthogonal  components 
at  the  appropriate  component  velocity. 

Since  constitutes  a  vertical  dipole  over  a  con¬ 
ducting  half  space,  the  fields  from  this  dipole  can  be  calcu¬ 
lated  using  the  Sommerfeld  integral  technique  (Wait  1970, 
Banos,  1966).  The  Hertz  potential  in  frequency  domain  due 

to  a  dipole  of  length  is  given  by 

z 


(4.11) 


where  ^2  “  +  (z-h)^  ^  from  Fig  4.7,  f(u)  is  the  Fourier 
transform  of  the  current  in  the  channel  (Equation  2.3)  given 
by 

and  the  value  of  k2  is  given  by 


^2  ~ 


(4/13) 


Using  the  Sommerfeld  integral  form,  equation  (4.11) 
can  be  written  as 
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where 


kjUi 


2 

1 


-0,{2+h)  . 

e  ^  J^(Xp)  dX 

A 


(4.15) 


,  m  =  1,2  (4.16) 

Sj 

kj  “[‘*>^^1^1  “  3  (4.17) 


and  J^(Xp)  is  Bessel  function  of  the  first  kind.  The  term 
2P  in  equation  (4.14)  is  a  correction  used  to  make  the  earth 
a  perfectly  conducting  surface. 

Wait  (1970)  approximated  the  value  of  the  integral 
given  by  equation  (4.15)  and  obtained 

P  =  (np)**  e”^^  erfc(ja)*'^)  -  (4.18) 

^1 


where 


P  = 


"’3^2^! 


and  0)  = 


A  is  given  by  the  relation 


“  Hq  (a^+joie^)  (4.19) 

Here  is  the  free  space  impedance  and  represents 
the  wave  impedance  of  the  ground  surface.  The  quantity  p 
defined  above  is  also  known  as  the  Sommerfeld  Numerical  Dis¬ 
tance. 

The  electric  field  can  be  calculated  from  the  Hertz 
vector  7(01)  =  [o ,  0 ,  by  the  equation 

^(03)  =  7y.¥(o))  +  7(01)  (4.20) 


Before  computing  the  field,  let  R  =  jr  -  r  |  where 
R  is  the  distance  from  the  observer  to  the  source  (R  is  either 
R^  or  Rg  in  Fig  4,7).  If  R^  =  17^1  represents  the  distance 
of  observer  from  the  center  of  the  source  and  also  the  dipole 
length  is  much  smaller  than  R^,  i.e. 


<  1 


then  R  can  be  written  to  first  approximation  as, 


R  =  R-+  Ir  -  r|  cos(|)  (4.21) 

Assumed  that  source  filament  is  electrically 
small  or  k2  <  <1 
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Now  to  find  the  field  due  to  a  finite  length  segment 
a  source  integration  technique  given  by  LeVine  and  Meneghini 


(1978)  is  followed.  The  field  components  from  a  vertical 
dipole  are  found  as 

( 


E  (01)  = 


joif  (aj)L  no 
2 


471k, 


sine  I  ’ 


-  (''"e'z)] 


.2 

3 _  e _ 

3r3z  ^2 


IP 


and 
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+  Sine  [ 
-JkjR^ 


R, 


-  2P 


E^(o))  = 


joif  (0))L2liO 


(4.22) 


47ik 


(&  •  >=“) 


2 

-jk2R2 


[sinc^Hk2L^(n2-(VR^)^)^ 


+  sinc^Hk2L^(n2-(VR^, 

•  “)] 


(4.23) 


R  '  is  the  same  as  R  but  for  the  filament  image. 

w  ^ 

T]„  is  the  ratio  of  speed  of  light  to  the  z  component  of  the 
z 
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and 


jcjf(u))L  uo 

E  (to)  - 5 —  Cos4> 


r^2 

>*2  / 


sine  -  ('?Re>x>] 


-  (TO’)^) 


(4.27) 


where  (J)  is  the  angle  made  by  r‘^  with  the  x-axis  and  the 
quantity  U  is  defined  as 


^  ^  2[i.jk  (z-^h)]  ^  r/^l 

r  3r  I  r  J 

The  value  of  V  can  be  found  from  the  equation 

2P  =  2  -  V 


(4.28) 


(4.29) 


where  P  has  already  been  defined  in  equation  (4.18), 

The  field  components  due  to  line  segment  can  be 
obtained  from  equations  (4.25)  to  (4.27)  by  rotating  the 
observer  location  parallel  to  the  ground.  Thus,  when  the 
field  components  are  known  then  their  vector  sum  provides  the 
electric  field  at  a  point  R  in  space. 

The  second  part  of  this  model  is  about  the  calculation 
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of  reflected  fields  from  ionosphere.  Only  single  bounce 
approximation  is  considered  because  the  subsequent  bounces 
contribute  very  little  to  the  total  field.  Fig  4.8  illustrates 
the  rays  contributing  to  the  single  bounce  reflection  model. 

Wait  (1962b)  has  calculated  reflection  coefficients 
for  incident  wave  on  ionospheric  layer  using  quasi  longtitudi- 
nal  approximation.  Ground  reflection  coefficients  are  also 
described  there.  According  to  these  calculations  and  as  shown 
in  Fig  4.9  and  4.10  the  ionosphere  effect  is  only  applicable 
up  to  100  KHZ  and  the  anisotropy  effect  begins  near  800  HZ  and 
is  maximum  near  100  KHZ. 

To  estimate  the  sky  wave  fields  simple  optical  ground 
reflection  technique  is  used  for  fields  due  to  each  line 
segment.  In  order  to  determine  the  extra  travel  time  and 
attenuation  due  to  increased  distance,  a  new  separation  vector 
from  the  segment  midpoint  to  a  new  observation  point  Rj  is  de¬ 
fined  as 


=  (Rx'  ^y'  2D)  -  K 


(4.30) 


where  D  is  the  height  of  the  ionosphere  above  the  ground. 

The  transmitted  field  Eq  measured  at  the  new  observa¬ 
tion  point  due  to  a  vertical  segment  is  given  as 


E°(ffi)  -  r  exp(jkRj)/R 


Frequency  vs  Reflection  Coefficient  for  Preserved  Polarization 
When  Waves  Reflected  from  Ionosphere  (Gardner,  1980) 


10*  10^ 
FREQUENCY  (HERTZ) 


O* 


Fig  4.10  Cross  Polarization  Reflection  Co-Efficient  from 
Ionosphere  as  a  Function  of  Frequency 
(Gardner,  1980) 
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si„c(-2^ 


(4.31) 


where  1  and  1'  are  unit  vectors  along  the  filament  and  its 

image,  respectively.  denotes  the  reflection  coefficient 

S 

due  to  the  ground.  Quantitites  having  superscript  (')  are 

due  to  the  image.  In  equation  (4.31)  the  sine  function  is  due 

to  the  phase  change  along  tbfe  filament.  From  the  transmitted 
— D  — s 

field  E  ,  the  sky  wave  field  E  's  on  the  g-round  are  given  by 
the  equations. 


(4.33) 


fan 


E®(u))  =  jlRlI  E°(aj)  +  llRlI  E^(a))  cossj 

^1+Rgllj  (4.34) 

Here  9  is  the  angle  between  the  wave  and  the  vertical 
and  symbol  |  jR  indicates  the  reflection  coefficient  for  the 
incident  wave  with  polarization  parallel  to  the  plane  of  in¬ 
cidence  to  the  reflected  wave  with  polarization  perpendicular 
to  the  plane  of  incidence. 

From  the  foregoing  discussion  and  Fig  4.9^  it  is 
obsious  that  radiations  up  to  1000  HZ  can  be  treated  as  pro¬ 
pagation  i'l  the  earth  ionosphere  waveguide.  The  geometry 
of  such  a  waveguide  configuration  is  shown  in  Fig  4,11.  Wait 
(1970)  developed  this  model  using  three  layer  Sommerfeld 
Integrals.  Using  vertical  line  sigment,  this  model  gives 

QSOS 


where 


M(a)) 


^(w) 


and 
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Fig  4.11  Geometry  for  Waveguide  Model  (Gardner,  1980) 
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are  the  reflection  coefficients  from  the  ionosphere  layer 
due  to  parallel  polarization  and  c^'s  are  the  solutions  to 

Rg(c)  Rj^(c)  exp  (jkgCD)  =  1  (4.36) 

®o^^^  is  Hankel  function. 

Equation  (4.35)  provides  relationship  between  the 
ionosphere  reflection  coefficients  and  the  ground  reflection 
coefficients.  Quantity  fn(z)  mentioned  in  equation  (4.35)  is 
defined  as 


exp(jk2C  z)  +  Rg.(c  )  exp(-jk2C  z) 

f  (z)  .  - 2-2 - i-,-2 - ±JS—  (4.37) 

[  2  f 


From  the  Hertz  potential  of  equation  (4.35)  the  elec 
trical  field  due  to  a  vertical  filament  is  calculated  using 
equation  (4.20),  Therefore 


E 

z 


(0))  =  ( 


(k2SnP)f„(h)f„(z)5„Cc) 
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(4.38) 


Also  the  field  from  the  horizontal  filament  can  be 


derived  from  equation  (4.38)  by  using  the  reciprocity  theorem 

(Wait  1960).  Since  the  fields  due  to  E  (co)  and  E  (u)  are 

X  y 

negligible  when  the  observation  point  is  near  ground,  the 

net  field  due  to  E_(a))  is  given  by 

z 


r 

L  j4ue^h 


6 


n 


(4.38) 


Here  f„  and  c  are  the  same  as  defined  by 
n  n 

(4.36)  and  (4.37).  However,  is  the  solution 


equations 

to 


*„(!>) 


(4.39) 


Zq  mean  the  value  of  z  in  free  space.  In  fact,  the 
suffix  0  wherever  it  appears  refers  to  free  space  quantities. 
(See  Fig  4.11).  From  equations  (4.31),  (4.35)  and  (4.38)  it 
is  apparent  that  the  frequency  spectrum,  and  thus  the  value  of 
radiated  field,  depends  upon  (a)  finite  size  of  the  channel, 
(b)  velocity  of  propagation  of  current  pulse,  (c)  channel 
tortuosity  and  (d)  ground  conductivity.  A  lightning  return 
stroke  model  which  is  not  based  upon  these  four  important 
parameters  will  not  predict  EM  fields  accurately.  The  iono¬ 
spheric  reflection  discussed  in  this  model  is  only  important 
for  times  exceeding  500  usec  or  so,  otherwise  this  parameter 
may  be  neglected. 


The  results  of  this  model  are  compared  to  Lin's 
(1978)  data  in  Fig  4.12. 

Strawe  et  al.  Model  (1979) 

This  is  a  nonlinear  transmission  line  model  which 
accounts  for  consistent  return  stroke  current  and  channel 
resistance  as  a  function  of  time.  The  model  is  based  on  the 
idea  that  the  return  stroke  channel  diameter  increases  by 
several  orders  of  magnitude  in  about  the  first  100  microseconds. 
Thus,  due  to  this  time  varying  parameter,  the  other  parameters 
during  a  return  stroke  also  vary  with  time.  Particularly  the 
channel  resistance  drops  to  a  great  extent  and  as  a  conse¬ 
quence  the  current  increases  in  the  same  order.  In  fact, 
not  only  the  channel  diameter  changes  with  time  but  para¬ 
meters  like  pressure,  temperature  and  particle  density  also 
change  with  time.  All  these  parameters  interact  during  the 
return  stroke  process  and  affect  the  size  ^nd  shape  of  the 
current  pulse. 

The  current  calculations  in  this.model  are  based  on 
time  varying  resistance  per  unit  length.  Although  inductance 
and  capacitance  per  unit  length  also  vary  with  time  due  to 
corresponding  changing  channel  parameters,  such  changes  are 
considered  weak  functions  of  time  and  are  not  considered 
(Strawe  1979,  80),  However,  capacitance  variance  is  com¬ 
pensated  by  considering  large  channel  radius  instead  of  sub¬ 
centimeter  current  carrying  core.  The  channel  tortuosity 
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and  corona  sheath  charge  storage  effects  are  accomplished  by 
modification  to  the  hot  channel  velocity  factor  (vj), 
characteristic  impedance  (Zq)  and  by  increasing  the  resis¬ 
tance  per  unit  length,  R  (Strawe  et  al.  1980). 

The  model  representation  is  divided  in  two  different 
forms.  The  first  form  consists  of  uniform  transmission  line 
segments  as  shown  in  Fig  4.13.  Here  the  channel  height  is  2 
kilometers  (Strawe  1979).  The  switch  shown  in  the  channel 
closes  when  the  stepped  leader  merges  with  the  upward  posi¬ 
tive  leader  from  the  ground.  The  channel  resistance  is 
measured  for  each  segment  using  Braginsikil's  (1958)  theory 
of  the  development  of  a  spark  channel.  This  process  involves 
the  past  current  time  history. 

According  to  Braginsikil's  theory  the  return  stroke 

channel  carries  the  current  in  a  narrow  path  formed  of  a 

2 

high  temperature  gas.  As  the  Joule  heat  (i  R)  is  released, 
the  pressure  in  the  channel  increases  and  so  does  the  channel 
radius.  This  inflated  part  of  the  channel  in  turn  acts  as  a 
piston  on  the  remaining  gas  which  expands  with  supersonic 
speed.  This  gas  then  causes  a  shock  wave  in  the  channel.  The 
temperature  of  the  gas  near  the  shock  wave  is  much  higher 
than  the  temperature  of  the  still  gas,  whereas  below  the 
shock  wave  the  temperature  is  many  times  higher.  This 
phenomena  causes  low  gas  density  in  the  channel  and  defines 
the  piston  boundaries.  With  increased  temperature  and  re¬ 
sulting  ionization  in  the  channel,  the  electrical  conduc¬ 
tivity,  which  is  proportional  to  T^/^,  increases  (Braginsikii 
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b) RETURN  STROKE  GEOMETRY 


Fig  4.13  Stroke  Channel  as  a  Transmission  Line 

and  Its  Geometry  m  Strawe's  Model.  (Strawe,  1978) 


1958) 


The  channel  resistance  per  unit  length  is  measured 
from  its  conductivity  and  the  arc  radius.  For  time  dependent 
conductivity  a(t)  and  arc  radius  a(t)  the  channel  resistance 
per  unit  length  R(t)  can  be  obtained  from  the  following 
equations 


a^(t) 


(4.40) 


_ ^ 

R(t) 


=  -T  a(t) 


a^(t) 


(4.41) 


V 


In  equation  (4.40)  the  quantity  a^  defines  the  ini¬ 
tial  condition  of  the  current  path  radius.  Equation  (4.40) 
is  a  solution  of  Braginsikii's  arc  model,  which  will  be  pre¬ 
sented  in  the  following  paragraphs.  Equation  (4.41)  is  a 
form  of  ohm's  law. 

In  arc  model,  the  channel  geometry  is  assumed  lo¬ 
cally  of  cylindrical  symmetry.  The  energy  deposit  in  the 

2 

channel  is  given  by  current  local  time  history  (i  R).  Typi¬ 
cal  arc  channel  radial  profile  is  shown  in  Fig  (4.14a).  The 
idealized  relation  between  temperature  T  and  arc  radius 
assumed  by  Braginsikii is  shown  in  Fig  (4.14b).  In  fact, 
Braginsikii assumes  three  radial  regions  in  the  channel.  The 
central  "channel"  region  is  the  hottest  region  with  low  den¬ 
sity  plasma.  The  entire  shock  system  internal  energy  is 
contained  here  except  for  the  kinetic  energy  of  translation. 
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a)  TYPICAL  ARC  CHANNEL  RADIAL  PROFILES 


b)  BRAdNSKil'S  MODEL  IDEALIZATION 
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Fig  4.14  Arc  Channel  Radial  Geometry  (Strawe,  1979) 
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Outside  this  region  is  the  "shell"  region  where  most  of  the 
gas  expelled  from  the  main  arc  is  trapped.  Although  this 
region  is  of  negligible  thickness,  almost  all  the  kinetic 
energy  of  the  shock  is  contained  here.  The  outermost  region 
is  known  as  "ambient  region"  and  it  is  beyond  the  shock  or 
shell  radius  a(t).  This  region  does  not  contribute  to  or 
take  energy  from  the  channel.  In  fact,  all  the  energy  is 
conserved  by  the  shock  and  shell  region  as  assumed  by 
Braginsikii.  The  relative  positions  of  these  regions  are  shown 
in  Fig  (4.15). 

From  piston  approximation,  as  mentioned  earlier,  if 
Oq  is  the  ambient  air  mass  density,  then  in  general,  the 
channel  pressure  p  is  calculated  from 

(P-Pq)  =  kp  pQ  ^  (4.42) 

where  p^  is  the  ambient  pressure  and  kp  is^the  efficiency 
factor  (kp<l).  Now,  if  p  is  much  larger  than  p^  then 
equation  (4.42)  reduces  to 

p  =  kp  Pq  (4.43) 

The  above  equation  provides  pressure-velocity  re¬ 
lation^  and  the  momentum  of  the  system  is  conserved. 

To  describe  the  channel  development,  the  energy  con¬ 
servation  equation,  the  thermal  and  caloric  equations  are 
formed.  The  energy  balance  equations  are 


Fig  4.15  Radial  Regions  in  Braginsikii ' s  Arc  Channel 
Model 
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where  W  M  and  €  are  internal  energy,  channel  mass  (Tra  p) 
and  internal  energy  per  unit  mass,  respectively.  and 
represent  rates  of  heat  conduction  and  heat  radiation  from 
the  lightning  channel  to  the  outer  shell  (Fig  4.15).  Equa¬ 
tion  (4.45)  elucidates  the  phenomena  of  mass  transfer  with 
heating  of  the  channel. 

Thermal  and  caloric  equations  from  thermal  physics 

are 


=  (Z-H)PkT 


(4.46) 


e  =  3£  +  =  - R— 

P  m^y  (Y'l)P  (4.47) 

Here  Z  is  the  average  ionic  charge,  whereas  p,  k,  T, 
and  m^^  are  channel  mass  density,  Boptzmann's  constant, 
channel  temperature  in  degree  kelvins  and  average  atomic 
mass  (kg),  respectively.  Here  y  is  the  effective  ratio  of 
specific  heats.  BraginsikfL  assumes  second  form  of  equation 
(4.47)  with  constant  value  of  y  He  also  assumes  that  a 
is  independent  of  time  and  temperature  and  is  a  constant 
quantity  (assumed  a  =  2.22x10*^  mhos/meter).  Thus,  under 
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these  assumptions  the  channel  resistance  per  unit  length 
becomes  the  function  of  a(t). 


Since 


W^=eM  =  (TTb^a") 


(4.48) 


Substituting  equations  (4.43)  and  (4.47)  gives 


*2 

C  (Y-  1) 


(4.49) 


After  another  substitution  of  equation  (4.44)  with 
equations  (4.41),  (4.43)  and  (4.49),  the  result  becomes 
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2,  r  d  2.2  ^  .2  d  ,2.] 
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^  P  °  ^  2(Y-2)a2J 


(4.51) 
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Let  E  =  1  + 


a2 


2(Y-1)&' 


(4.52) 


Apparently  E  is  a  function  of  time.  But  it  is  assum¬ 
ed  by  Braginskii  that  E  is  constant  due  to  power  law  and  hence 
equation  (4.51)  becomes 
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Integration  of  equation  (4.54)  gives 


a2  ,  -  2 


■'  i!)  1/3 


where 


(4.53) 


(4.54) 


(4.55) 


8  =  V2Ti2k  p^E 
p  o 


Obsiously  equation  (4.40)  and  equation  (4.55)  are 
identical . 

Since  the  channel  radius  is  a  funciton  of  time 
(assumed  a  simple  power  form)  therefore 


a 


(4.56) 


L 


75 


Also  channel  current  depends  on  resistance  per  unit  length 
and  some  function  of  time, therefore  'i'  can  be  expressed  as 


i  =  Rt‘ 


(4.57) 


From  equations  (4.52)  and  (4.56)  the  result  after  simplifica¬ 
tion  becomes 


E  =  1  + 


(Y  -  1) 


(4.58) 


Substitution  of  the  values  from  equations  (4.57)  and 
(4.58)  in  equation  (4.53)  after  comparing  the  powers  of  t 
results  in 


2n+3 


(4.59) 


Braginskii  assumed  that  the  valu^  of  a^  is  very 
small  as  compared  to  a(t).  Substitution  of  equation  (4.57) 
in  equation  (4. 55)  yields 


2 


The  solution  of  this  integral  gives  a  result  of  the 


form 


,1/3  ^1/2 


(4.60) 


^6 


where 


"  (2n+3)^'^ 

This  is  Braginskii  analytic  solution  for  the  arc 
model  of  lightning  channel.  However,  as  suggested  by  Strawe 
(1979)  the  numerical  solution  of  equation  (4.40)  is  convenient 
and  he  has  used  this  equation  in  his  first  model. 

The  second  model  suggested  by  Strawe  (1979)  is  a 
refined  version  of  the  first  model.  This  last  model  accounts 
for  the  effective*  energy.  Also  a  is  no  longer  considered  as 
constant.  Instead  values  of  a  are  computed  as  a  function  of 
channel  temperature  and  pressure.  Equation  (4.42)  is  re¬ 
written  by  applying  Newton's  Second  law  to  the  shell  with 
ambient  and  shell  pressures  working  in  opposite  direction. 

Thus 


(p-p  )  2it  ^  “  Time  rate  of  change  of  momentum 


^  ^  shell  mass  x  aj 


(4.61) 


Since  the  shell  does  not  exist  before  the  arrival  of 
the  shock  wave,  therefore,  initial  mass  of  the  shell  is  con¬ 
sidered  to  be  zero.  Equation  (4.45)  provides  a  relation 
between  energy  and  mass  transfer.  Q^.  is  considered  as 
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as  ultraviolet  radiation  from  the  channel  which  is  absorbed 


by  the  shell  to  produce  plasma  in  the  channel.  Evidently, 
all  the  ultraviolet  radiation  is  not  absorbed,  and  some  of 
it  is  wasted.  Taking  this  into  account,  Q^.  can  be  repre¬ 
sented  as 


Sr  ^total  ^a 


where  is  the  total  radiation  and  f^  is  the  absorp¬ 

tion  factor.  Therefore,  equation  (4.44)  can  be  modified  as; 


IT  a^(t) 


^total^^'^a^ 


(4.62) 


It  is  assumed  in  this  model  that  the  total  radiated 
power  is  the  transparent  radiation  and  proportional  to  the 
radiation  factor  and  channel  mass.  Thus, 

where  Fj,(t)  is  the  radiation  factor  of  the  arc  channel. 

The  heat  conduction  term  has  the  value  given  by 

=  2tt  mK^  (4.64) 

where  m  has  value  2  or  3  and  it  changes  with  a(t).  K.p  is 
a  function  of  T  and  p  and  it  is  related  by  the  Weedmann- 
franz  law  for  values  of  T  exceeding  15,000  K  (F.  Reif,  65). 


The  values  of  a(p,T)  and  K,p  used  in  second  model  are 
a  curve  fits  to  the  results  calculated  by  Plooster  (1971) 


_  1.52x10^ 

-  -  13.4  ^  .-54 


(4.65 


Kt  (P.T)  = 


\  (oTkIO""^)^ 

10  29.1 


3  4 

1  +  10^  (1~) 


(4.66 


The  current  waveform  used  in  this  model  is  similar 
to  Bruce-Golde  model  waveform  with  peak  value  of  20  kA  at 
1.5  usee  and  a  half  value  fall  time  of  40  psec.  The  result¬ 
ing  temperature,  pressure  and  radius  of  the  channel  due  to 
this  current  are  shown  in  Fig  4.16. 

Em  fields  from  the  current  are  calculated  in  time 
domain  with  the  assumption  that  earth  is  plane  and  perfect 
conductor . 

The  comparison  of  models  presented  in  Chapter  III, 


.V. 


and  IV  will  be  discussed  in  Chapter  VIII. 


V.  Electromagnetic  Fields  Due  to  Verticle 
Channel  of  Variable  Height  and  Having 
Arbitrary  Current  Flow 

Most  of  the  models  presented  in  the  last  two  chapters 
predicted  EM  fields  at  ground  level  due  to  a  vertical  channel 
having  arbitrary  current  flow  i(z',t).  Master  et  al.  (1981) 
have  introduced  a  lightning  return  stroke  model  in  which  the 
return  stroke  electric  and  magnetic  fielu  measured  above  the 
ground  level  can  be  reproduced.  The  equations  which  predict 
EM  fields  in  space  have  not  been  derived  in  their  presentation 
(Master  et  al.  (1981),  However  in  this  chapter  we  will  derive 
these  equations  of  electromagnetic  fields  to  have  an  analytic 
insight  of  the  model  and  it  is  hoped  that  these  derivations 
will  enable  the  reader  to  understand  better  the  modifications 
made  by  the  author  in  the  original  Master  et  al.  model.  The 
model  and  the  modifications  will  be  discussed  in  the  next 
chapter. 

The  geometry  of  a  straight  vertical  channel  of  height 
H  above  a  perfect  conducting  ground  plane  is  shown  in  Fig  5.1. 
In  the  figure  the  observation  point  P(r,(j),z)  is  at  attitude  z 
and  a  vertical  dipole  of  length  dz'  is  at  height  z'.  Equa¬ 
tions  for  EM  fields  will  be  derived  using  the  method  similar 
to  Uman  et  al.  (1975).  From  electromagnetic  field  theory, 


we  know  that 


wZVi' 


and 


7x(E+||)  =  0 

where  S’  is  the  vector  potential,  Equation  (5,2)  holds  when 

E  +  =  -V(l)  (5.3) 

Equations  (5.1)  and  (5.3)  can  predict  time  dependent 
fields  provided  value  of  S  is  known  for  the  point  in  space. 

The  Lorentz  condition 


VA  = 


-n 


o 


(5.4) 


in  conjunction  with  Maxwell's  equations  lead  us  to  a  result 

-  “o  (®-5) 

The  solution  of  the  inhomogeneous  wave  equation  (5.5) 
is  very  similar  to  the  solution  of  Poisson's  equation  (Plonsey 
and  Collin,  1961).  Note  that  the  resultant  value  of  vector 
potential  in  free  space  is  only  dependent  on  the  current  den¬ 
sity  T.  Thus,  if  current  information  in  the  lightning 
channel  is  known,  then  the  value  of  S  is 
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where  C  =  l/Zy^  and  dV  is  infinitesimal  volume  having 
current  density  T.  Fig  5.2  illustrates  the  geometry  in  the 
computation  of  retarded  potentials  given  by  the  equation  (5.6) 
To  implement  these  results  in  the  case  of  a  vertical 
channel  as  shown  in  Fig  5.1  and  meet  the  boundary  conditions 
an  imaginary  current  of  equal  amount  must  exist  vertically 
below  the  channel. 

Now  consider  a  dipole  of  length  dz*  along  the  verti¬ 
cal  channel  at  height  z’  (Fig  5.1).  Using  equation  (5.6), 
the  differential  vector  potential  at  point  P  due  to  dz‘  is 

dS(r,(D,z)  =  i  (z’,t-R/c)  dz’  (5.7) 

In  spherical  coordinates 

a^  =  a^^cose  -  ag  sin0 
Hence,  equation  (5.7)  becomes 

die  *  ^  i(z',t-R/c)  ^  a^cosG-agSine  J  dz'  (5.8) 


But  from  equation  (5.1) 


Fig  5.2  Geometry  for  the  Calculation  of  Retarded  Potential 
at  Field  Point  Due  to  Charge  Location  in  Volume  v* 
(Uman  et  al..  1969) 


dn  =  7xdlf 


(5.9) 


In  polar  coordinates  the  curl  of  X  is 

[  ilFe*  “ 

Equations  (5.8),  (5.9)  and  (5.10)  yield 


but 


8  j  x  I  ■.._R\  _  _1  3 

3^  i(z  .t-3)  -  -^ 


i(z',t-f) 


Therefore, 

4  i(z',t-|)]  dz'  (5.11) 

R 

Since  the  image  element  (-a^dz')  has  the  current  i(z',t) 

in  the  a  direction,  the  resultant  magnetic  field  d5  at  Point 
z 

P  is  dependent  on  vector  distance  (subscript  m  indicates 
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quantities  due  to  imaginary  element).  Using  the  result  from 
equation  (5.11),  we  get 


m 


uosine„  1 

1 

3i 

4Tr  1 

c5~" 

m 

• 

R  . 

^  . 

i(z' 

m 

*  c  ‘ 

-)  + 


(5.12 


where  is  the  angle  made  by  R  with  z-axis. 
m  °  "  m 

The  resultant  magnetic  field  dS  at  point  P  is  given 
by  the  vector  sum  of  equations  (5.11)  and  5.12),  therefore. 


‘^^total 


=  dS  +  dS 

m 


_  UO 
■  ?? 


It 


R 


L  ro 


fi .  dz ' 
<l> 


(5.i: 


Let  r  be  the  vertical  distance  of  point  P  from  the  channel 
(Fig  5.1).  Values  of  R,  R^,  Sine  and  Sine^^^  are  given  by 


R  = 


-F 


^+(z-z’ )^ 


(5.14; 


Sine 

r^+(z-z’)2 

(5.16) 

Sin9„ 

m 

r^+(z+z')^  j 

(5.17) 

where  z  is  the  vertical  height  of  the  point  P, 

Substitution  of  these  values  in  equation  (5.13)  yields 


The  total  magnetic  field  due  to  .the  channel  current 
length  l(t),  which  is  a  function  of  time  and  is  given  by 
equation  (3.9),  is  found  by  integrating  equation  (5.27)  from 
zero  to  i(T).  Therefore, 
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To  find  the  electric  field  due  to  an  infinitesimal  cur¬ 
rent  dipole  dz'  at  a  height  z’  along  the  channel,  equations 
(5.3)  and  (5.4)  are  used.  From  equation  (5.3) 


dE  +  Iy  dl  =  -  ?d<t 


(5.20) 


and  from  equation  (5.4) 


d  $  =  -c' 


V.  dl 


dt 


(5.21) 


In  polar  coordinates  the  divergency  of  A  is 


V.  I  ■  Ir  (A.) .  1 


2  3r  ''  rSinO 


§  [I?  (*9  ^  (5.22) 


Determining  the  divergence  of  equation  (5.8)  yields 


jX  _U0C0S9  fS  //O  S  r  r,  t  •  _ 


- ±(z*  ,t-^)  fr(Sin‘e)dz' 

4tvR'^  Sine  ^ 


uocose 

4itR^ 


Cos9  i(z'/t~)dz' 
c 


uocose 

4ttRc 


i(z,t~)  dz' 


(5 


Substitution  of  equation  (5.23)  in  equation  (5.21)  yields 

t 


d<Ii  = 


cR 


i(z 


dz ' 


(5 


because  c  =  l/y^e^. 


The  negative  gradient  of  equation  (5.24)  gives  the 


result 


-  7d$  = 


1  3 

i  (z ' 

t 

sine  r 

1 

1  .f  ^  ff  _ 

4iTeQR  [ 

? 

/  i(z  ,T 

t-s; 

c ' 

1  ^ei'*^' 

(5 


To  find  ^  dK,  take  the  derivative  of  equation  (5.8).  The  re¬ 


sult  is 


=  il£  f 

4Tr 


Sin0  3 
R  ft 


(5.26) 


Substitution  of  equations  (5.25)  and  (5.26)  in  equa¬ 


tion  (5.20)  gives 


=  4ife  [ fa  / 

Slnsf^y  i(z',T-|)dT  +^l(z',t-|) 


*  ^  h  ■'-!)]  *e  . 


(5.27) 


The  field  due  to  an  infinitesimal  imaginary  element 

at  a  distance  from  the  observation  point  P  (Fig  5.1)  can  be 

found  from  equation  (5.27)  by  replacing  R  with  R  .  6  with  e  . 

ra'  m' 

“'“h  \ 

Therefore,  the  electric  field  dE  is 

m 
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+  Sine 


i(z 


+ 


1 


i(z 


(5.28) 


) 

Before  finding  the  total  electric  field  by  adding 
equations  (5.27)  and  (5.28)  it  is  convenient  to  convert  the 
differential  fields  to  cylindrical  coordinates.  From  Fig  5.1 
it  can  be  found  that 


=  Sine  dj.  +  Cose  d^  (5.29) 

d  =  Cose  -  Sine  (5.30) 

=  cos0^  +  cos0„  (5.31) 

K_  m  r  m  z  '  ' 

m 

a  „=  Cos0„  a^  -  Sine^  &  (5.32) 

m  m  r  m  z  ^  ' 


Substitution  of  these  values  in  equations  (5.27)  and 
(5.28)  results  in 


dE  = 


dz ' 


(.acosVs.i^ei  j 


*  RC°.s.V3.1«^ 


SitA 

C^R 
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.  3Cos0Sin0  w_i  *  R 
+  - K -  i(z  ,  t  -  - 

cR^  ® 


Cos0Sin0  3  ^  ,  *  rJ  . 

-  It 

4 


(5.3 


A  similar  equation  for  dE^  can  be  determined  ex¬ 
cept  that  0  and  R  are  replaced  by  9^^^  and  R^j,,  respectively. 
Similarly  to  equations  (5.16)  and  (5.17)  the  values  of  Cos 

and  Cos9„  are  defined  as 
m 


Cos0  = 


2  -  Z 


(5.3 


Cos0 


2  +  Z  ' 


m  R_ 


(5.3 


Using  equations  (5.16),  (5.17),  (5.34)  and  (5.35), 
the  resultant  field  at  point  P  is 


=  dE  +  dE„ 
total  *  m 


2(z-zM^-r^ 

»5 


i(z',T  -  g)dT 


CR^  ° 


L'***;  V- 


i 


m  m 

Jmipr  r  .3(z-z^)r  t(^.,t.|) 

V  R  1  cR 

^  ^  It  ^  ^  J 


c'  p5 


+  2i|!£l2l  1, 


m  ■* 


(5.36) 


The  total  electric  field  due  to  channel  height  lor  2 
is  obtained  by  integrating  equation  (5.36)  and  the  result  is 


^  4 


l(t)  > 

dE 


(5.37) 


k 


N 


Values  of  R  and  R^j^  in  equation  (5.36)  have  already  been 
defined  in  equations  (5.14)  and  (5.15).  Equation  (5.36)  is 
exactly  the  same  as  Master  et  al.  (1981)  have  used  in  their 
model  without  deriving  it. 

Note  that  the  equation  (5.36)  depicts  the  total  field 
comprising  of  radial  and  vertical  components.  However,  the 
airborne  measurements  of  the  electric  field  is  only  in  verti- 

A 

cal  direction;  whereas  component  lies  in  the  plane  of 

the  steady  state  aircraft,  and  its  value  is  not  recorded. 

Equation  (5.19)  is  the  predicted  value  of  magnetic 

field  in  d,  direction.  In  experimental  data  the  value  of  F 
<P 

is  recorded  in  two  orthogonal  components.  These  components 
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are  along  the  nose-to-tail  and  wing-to-wing  direction  as  shown 
in  Fig  5.3.  These  components,  along  with  the  aircraft  bearing, 
are  used  to  find  the  relative  direction  of  the  lightning  re¬ 
turn  stroke.  The  detailed  discussion  on  the  validity  of  this 
data  and  direction  finding  will  be  discussed  in  one  of  the 
next  chapters.  However,  it  is  apparent  from  equation  (5.19) 
that  in  ideal  conditions  correct  bearing  of  the  lightning 
channel  with  respect  to  aircraft  can  be  found. 


VI.  Master  et  al.  Model  (1981) 
And  Modifications 


Master  et  al.  Model  (1980) 

This  model  is  a  modified  version  of  the  Lin  et  al. 
model  (1980)  and  provides  reasonable  estimates  of  the  return 
stroke  fields  in  space.  As  suggested  by  the  authors  of  this 
model,  the  EM  fields  computation  was  restricted  to  an  alti¬ 
tude  of  10  km  and  at  ranges  from  20  m  to  10  km.  However, 
Master  et  al.  did  not  specify  explicitly  the  reasons  for 
restricting  these  ranges. 

Similar  to  the  Lin  et  al.  model,  the  return  stroke 
channel  is  like  a  vertical  thin  wire  antenna  of  height  H  above 
a  perfectly  conducting  ground.  The  geometry  of  the  model  is 
shown  in  Fig  5.1  and  the  mathematical  expression  that  predicts 
the  fields  above  ground  level  were  derived  in  Chapter  V  and 
were  given  by  equations  (5.19)  and  (5.37).  The  current 
parameters  of  this  model  are  due  to  subsequent  return  stroke 
rather  than  to  a  first  return  stroke.  Subsequent  return 
strokes  were  considered  because  they  have  fewer  branches 
and  relatively  constant  upward  propagation  velocity. 

The  modifications  to  the  Lin  et  al.  model  (1980)  are 
based  on  defining  the  current  parameters.  The  Lin  et  al. 
model  assumes  that  a  short  duration  breakdown  pulse  current 
propagates  upwards  with  constant  amplitude  and  waveshape. 


But,  in  this  model  the  breakdown  pulse  is  assumed  to  decrease 
with  height  because  the  subsequent  stroke  initial  peak 
luminosity  changes  considerably  with  height  (Jordan  and 
Uman  1980).  The  attenuation  in  the  breakdown  pulse  also 
removes  the  mirror  image  problem  observed  in  the  Lin  et  al. 
model.  Mirror  image  is  the  field  change  of  opposite  polarity 
to  that  of  the  initial  field  when  the  breakdown  pulse  reaches 
the  top  of  the  channel  and  this  effect  is  almost  never  ob- 
served  in  the  subsequent  return  stroke  (Lin  et  al.  1980). 

The  channel  current  waveform  consists  of  three  com¬ 
ponents  called  the  uniform  current,  the  corona  current, 

i  and  the  breakdown  pulse  current,  i  .  Mathematically  it 
c  p 

can  be  represented  as 

l(z,t)  =  !„  +  ip  (8.1) 

A  graphical  comparison  of  the  currgnt  waveforms  in 
the  Lin  et  al.  model  (1980)  and  Master  et  al.  model  (1981) 
are  shown  in  Fig  (6.1).  It  is  obvious  from  Fig  6.1  that  the 
minimum  value  of  the  corona  and  breakdown  pulse  are  equal  to 
the  uniform  current  in  the  Lin  et  al.  model;  however,  this 
value  is  considered  as  zero  in  Master  et  al.  model  (1981). 
Although  not  shown  in  Fig  6.1,  the  rise-to-peak  of  the  break¬ 
down  pulse  has  also  been  modified  from  that  of  Lin  et  al. 
(1980)  so  as  to  make  these  parameters  consistent  with  the 
results  of  Weidman  and  Krider  (1978). 
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Pig  6-1  Graphical  Comparison  of  Return  Stroke  Channel 
Current  Waveforms  (a)  L  et  al.  Model  (1980) 
(b)  Master  et  al.  Model  (1981) 


Except  for  the  modifications  discussed  above,  the 
Master  et  al.  model  (1981)  uses  all  the  current  parameter 
described  in  Lin  et  al.  model  (1980)  in  Chapter  III.  In 


Master  et  al.  model  the  height  dependent  attenuation  of  the 

-2 '  /  X 

breakdown  pulse  current  takes  the  form  ip  e  '  where  X  = 
1500  meters.  Thus,  the  breakdown  pulse  decreases  with  height 
in  exactly  the  same  manner  as  does  the  corona  current  in  Lin 
et  al.  model  (1980). 

Modifications 

During  this  research  only  two  modifications  in 
Master  et  al.  model  were  studied.  These  modifications  are 
related  to  the  variable  propagation  velocity  and  the  wave¬ 
shape  of  the  corona  current.  Master  et  al.  model  (1981) 
is  in  fact  based  on  the  statistical  averages  of  various  para¬ 
meter  in  a  return  stroke.  But  the  parameter  chosen  for 
modification  are  such  that  their  statistical  average  will 
considerably  effect  the  estimated  value  of  the  EM  fields 
above  the  ground. 

Lin  et  al.  (1980)  have  mentioned  that  error  in  the 
measurement  of  propagation  velocity  could  be  as  high  as  30 
to  60  percent.  In  fact,  measured  velocities  depend  on  the 
height  above  the  ground,  at  which  the  measurement  was  made. 
Also,  Uman  and  McLain  (1970)  have  proved  that  for  distance 
0  >  50  km 


i(D.t)  =  B(t+i)  for  1  <  H  (6.2) 


where  l,  H  are  the  instantaneous  and  maximum  length  of  the 
return  stroke  channel. 

Equation  (6.2)  provides  a  simple  relation  between  the 
channel  current,  velocity  of  current  pulse  propagation  and 
the  distant  radiated  magnetic  field.  Thus,  if  larger  return 
stroke  velocity  is  assumed,  then  the  current  peak  value  and 
the  waveshape  is  altered.  These  arguments  justify  that 
variable  velocity  should  be  considered  in  modeling  the  current 
waveform. 

Almost  in  all  the  nonlinear  models  the  velocity  of 
propagation  is  varying  because  of  nonlinear  parameters  L  and 
C.  Braginskii  (1958)  has  also  considered  the  nonlinear  velocity 
in  an  arc  channel.  Lin  et  al.  (1980)  has  suggested  that  refine¬ 
ment  in  the  model  can  be  achieved  by  taking  into  consideration 
the  decrease  in  velocity  with  height.  Jordan  and  Uman  (1980) 
have  found  that  the  peak  luminosity  of  the  subsequent  strokes 
decreases  to  half  the  initial  peak  in  less  than  one  kilometer 
above  ground.  In  recent  studies  Uman  et  al.  (1982)  report 

Q 

that  average  velocity  for  subsequent  stroke  is  1.2x10  m/sec  and 
this  velocity  decreases  with  height. 

The  author  has  found  out  during  these  studies  that  if 
the  attenuation  factor  of  the  current  pulse  in  Master  et  al. 
model  is  not  changed  but  the  velocity  of  propagation  is  in¬ 
creased  then  the  mirror  image  effect  can  be  found. 

It  was  also  found  by  trial  and  error  that  the  velocity 

Q  O 

variation  between  1.5x10  m/sec  and  .5x10  m/sec  oculd  be 
achieved.  Since  the  change  in  velocity  is  considered  height 


dependent  and  there  was  no  experimental  data  available  with 
the  author  to  put  velocity  in  this  form,  v^  =  v^  for 

the  first  stroke,  therefore,  velocity  variation  was  consi¬ 
dered  linear  instead  of  exponential.  Modification  in  sub¬ 
sequent  stroke  is  not  studied  because  experimental  data 
for  subsequent  stroke  was  not  made  available.  The  modifi¬ 
cation  made  in  Master  et  al,  model  is 

v(z)  =  Vq(1-.33z/H)  (6.3) 

8 

where  v^  was  taken  as  1.5x10  m/sec  and  H  =  5  km.  This 
relation  is  satisfied  for  a  return  stroke  travel  time  of 
50  ysec,  which  is  the  same  in  the  Master  et  al.  model 

Q 

with  constant  velocity  v  =  1x10  m/sec.  Berger  et  al, 

(1975)  have  shown  that  average  value  of  subsequent  dura¬ 
tion  has  exceeded  32  usee  for  50%  of  cases.  If  these  values 
are  taken  into  account,  then  the  propagation  velocity  can 
still  be  increased.  The  results  of  this  modification  are 
shown  in  Fig  8.16  to  8,18, 

The  second  modification  made  in  tht  Master  et  al. 
model  is  in  corona  current  waveshape.  Bruce  Golde  (1941) 
suggested  the  double  exponential  decay  waveform  and  in  the 
later  years  most  of  the  researchers  were  basing  their 
models  on  this  waveform,  Uman  (1969)  and  Lin  (1978)  made 
significant  changes  in  the  Bruce  &  Golde  waveform.  Uman 
(1969)  suggested  attenuation  in  the  corona  current  waveform. 
Lin  (1978)  divided  the  channel  current  into  three  components 
as  given  by  equation  (6,1),  Recent  studies  in  current  wave¬ 
form  were  carried  out  by  Berger  et  al.  (1975),  Garbanatti 
(1982),  and  Eriksson  (1978),  The  typical  return  stroke 


current  measured  by  these  authors  is  shown  in  Fig  6,2. 
According  to  these  studies  for  the  first  and  subsequent 
return  strokes,  current  drops  to  50%  of  peak  value  in  al¬ 
most  30  psec.  Fig  6.2(a)  shows  a  large  variation  in  the 
current  waveform  for  the  first  return  stroke.  However, 
waveform  A  in  Fig  6.2(a)  indicates  a  fast  decay  in  the 
current  even  after  30  usee.  The  current  parameters  suggested 
by  Uman  et  al.  (1981)  for  the  first  return  stroke  of  Master 
et  al.  model  (1981)  produce  a  current  waveform  much  similar 
to  this  waveform.  But  the  EM  fields  for  first  stroke  pre¬ 
dicted  in  space  by  Master  et  al.  model  (1981)  and  shown 
in  Fig  8.13  to  8.15  are  different  from  those  of  empirical 
data  as  shown  in  Fig  8,3  to  8,12,  Because  during  this 
study  empirical  data  only  for  the  first  return  stroke  was 
made  available,  therefore,  a  modification  to  the  current 
waveform  for  the  first  return  stroke  is  suggested.  The 
modification  in  current  waveform  is  also  suggested  because 
the  current  parameters  in  first  stroke  changes  signifi¬ 
cantly  with  time.  In  Chapter  VIII,  it  will  be  proved 
that  combined  effect  of  modified  velocity  and  current 
waveform  improve  considerably  the  peak  value  of  predicted 
fields. 

For  the  current  waveform  modification  values  of 
a  and  $  are  varied  with  time  so  as  to  achieve  maximum 
value  of  ^  during  30  sec  to  50  sec  time  and  then  mini¬ 
mum  value  of  ^  in  the  subseuqent  50  sec.  It  is  assumed 
that  total  observation  time  is  100  sec.  Since  Master  et 


al,  model  current  waveform  for  first  return  stroke  con¬ 
forms  with  waveform  A  of  Fig  6.2(a)  up  to  30  ysec,  there¬ 
fore,  modification  in  the  waveform  was  incorporated  30 
psec  after  the  initiation  of  the  return  stroke.  The 
value  of  a  and  B  Increases  logrithmically  as 


a  =  1x10  X  In 


20x10 


+  .6 


+  >  30  psec  (6.4) 


3  =  3xl0°x  In 


20x10 


+  .6 


+  30  psec  (6,5) 


when  't'  is  in  seconds. 

The  quantity  .6  in  the  above  equations  is  added  to 

have  numerical  value  1  for  the  quantity  within  the  brackets 

<" . 

at  30  psec.  Between  30  to  30  psec  after  the  initiation  of 
return  stroke  the  values  of  a  and  3  are  increased  to  1,5 
times  of  its  initial  value  and  thus,  rapi^i  decay  in  the 
current  occurs  during  20  psec  period.  However,  in  the 
subsequent  50  psec  only  1,47  times  the  values  of  a  and 
3  increases  and  thus,  slow  rate  of  change  in  current  during 
this  period.  The  results  of  this  modification  are  shown 
in  Fig  8.19  to  Fig  8.21. 


VII.  Validity  of  Airborne  Data 


Measured  values  of  electromagnetic  fields  and  related 
parameters  in  lightning  environments  are  vital  for  the  testing 
of  return  stroke  models  and  to  ascertain  the  amount  of  electro¬ 
magnetic  coupling  to  the  aircraft.  The  airborne  measurements 
of  fields  have  certain  complexities  which  are  comparable  to  the 
ground  based  measurements.  Since  in  this  thesis  we  are  con¬ 
cerned  mainly  about  the  testing  of  the  models,  the  problems 
associated  with  the  acquired  data  which  are  used  to  justify 
the  model  will  have  to  be  discussed.  The  sources  of  error  in 
the  data  are  classified  in  four  categories,  namely, 

1.  Data  collection 

2.  Direction  finding 

3 .  Range  measurement 

4.  Channel  height  measurement 

This  classification  is  based  on  the  system  instru¬ 
mentation  and  calibration  used  by  the  Flight  Dynamics  Labora¬ 
tory  for  the  collection  of  the  airborne  data  during  the  years 
1979  -  1981,  and  the  data  made  available  to  the  author  for 
the  comparison  of  the  models  presented  in  this  thesis. 

Data  Collection  Problems 

The  most  serious  problem  in  interpreting  the  airborne 
data  is  due  to  the  conducting  outer  surface  of  the  aircraft  on 


which  the  field's  sensors  are  installed.  The  surface  current 


on  the  aircraft  caused  by  the  incident  electric  field  distort 
the  electric  field  in  the  vicinity  of  the  sensor  and  thus,  the 
data  recorded  does  not  truly  indicate  the  incident  electric 
field.  Additional  uncertainties  are  added  due  to  aircraft 
orientation  with  respect  to  the  return  stroke  channel  and  the 
polarization  of  the  incident  E-field. 

Perala  et  al.  (1982)  have  provided  an  analytic  treat¬ 
ment  to  this  problem.  Let  Ej^(r,t)  and  Hj^(r,t)  be  the  incident 
fields  on  the  aircraft  having  surface  area  S  and  assume  that 
initially  the  aircraft  has  no  surface  charges  and  surface 
currents.  Then  the  incident  field  causes  a  surface  current 
density  J_  and  charge  a.  Due  to  this  current  and  charge,  the 
scattering  fields  E.(r,t)  and  H„(r,t)  are  then  produced. 

w  S 

The  total  fields  E,j,  and  Ep  in  the  vicinity  of  the 
aircraft  are  given  by 


E^(r  ,t)  =  Eg(r  ,t)  +  E^(r  ,t) 


(7.1) 


H^(r  ,t)  =  EgCr  •,t)  +  !r.(r  ,t) 


(7.2) 


The  current  density  J  and  a  are  given  by  the  equa 

tions 


^  s  '  ^  ^  ~  flxH  q,  {  »  i  ) 


(7.3) 


and  a(  r  ,t)  =  efi.f^(rg,t) 


(7.4) 
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where  the  r_  vector  is  on  the  aircraft  surface  and  fi  is  the 
s 

normal  to  the  surface  at  that  point. 

It  is  obvious  from  equations  (7.1)  and  (7.2)  that 
the  true  value  of  the  incident  EM  fields  due  to  lightning  are 
difficult  to  measure  even  using  inverse-scattering  techniques 
(Baum  et  al.,  1982).  Even  with  the  same  sensor,  the  fields 
measured  at  another  aircraft  under  the  similar  conditons  may 
have  significant  differences. 

Another  problem  in  airborne  data  collection  is  re¬ 
sonant  effect  of  the  aircraft.  It  has  been  observed  that  at 
certain  frequencies  the  a/c  body  causes  resonance  in  the  EM 
fields.  In  fact,  this  is  a  multi-components  damped  sinusoidal 
behavior  of  the  aircraft  body  which  distorts  the  incident  field 
and  corrupts  the  measured  data. 

The  EM  fields,  which  are  of  interest  in  model  testing 
comprise  the  vertical  component  of  the  E  field  in  9._  direc- 
tion  and  has  two  orthogonal  components  of  H_,  the  magnetic  field 

•  -f 

The  components  of  H.  are  measured  with  crossed  magnetic  loops 
and  then  used  to  find  the  angle  of  the  lightning  channel  with 
respect  to  the  aircraft  bearing.  In  actual  flight  conditions 
the  instantaneous  banking  and  pitch  angles  of  the  aircraft 
reduces  the  perpendicular  area  of  the  sensing  coil.  The 
reduction  in  area  is  equal  to  cosine  of  the  respective  angles. 
Although  the  data  were  collected  in  level  flights,  angular 
errors  can  be  caused  due  to  air  turbulence,  etc.  Crossed 
magnetic  loops  were  installed  along  the  nose-to-tail  symmetry 
and  wing-to-wing  direction.  A  typical  sensor  location  is 


shown  in  Fig  S.3. 

For  electric  field  measurement,  the  flush  plate 
dipole  (FPD)  were  used.  The  disadvantage  of  this  sensor  is 
due  to  the  attenuation  caused  by  the  added  series  resistance 
which  increases  the  overall  input  impedance  of  the  system 
(Baum  et  al.,  1982).  Calibration  results  of  this  sensor  in 
the  actual  system  were  within  10%  error  from  the  theoretical 
value  and  a  minimum  sinusoidal  signal  of  50  yV  amplitude  could 
be  detected.  Thus,  the  electric  field  sensor  could  detect  a 
minimum  reading  of  less  than  20  V/m  which  corresponds  to  a 
range  of  over  30  km.  The  magnetic  field  sensor  consisted  of 
one  cylindrical  moebius  loop  and  the  system  response  had  error 
within  5%.  The  detailed  theory  about  sensors  and  their  per¬ 
formance  is  given  by  Rustan  et  al.,  1982. 

Direction  Finding  Problems 

Crossed  magnetic  loops  were  used  tp  find  the  relative 
bearing  of  the  lightning  channel  with  respect  to  the  aircraft. 
It  is  assumed  that  the  channel  is  vertical  and  its  magnetic 

A 

field  is  only  in  the  a^  direction  (in  cylindrical  coordinate). 
The  voltages  induced  in  the  two  magnetic  field  loops  are  pro¬ 
portional  to  sin0  in  one  loop  and  cos0  in  the  other  loop. 

The  relative  bearing  of  the  lightning  channel  is  shown  in 
Fig  7.1  and  the  two  loops  are  supposed  to  be  identical  in  shape 
and  size.  The  ratio  of  the  loop  voltages  gives  the  value  tan0 
and  thus,  the  value  of  0  defines  the  relative  bearing  of  the 
channel.  Since  in  actual  measurements  the  magnetic  field 


Fig  7.1 


Relative  Bearing  of  the  Li| 
with  Respect  to  Aircraft  A: 


components  were  used,  the  sign  of  the  magnetic  fields  should 
determine  the  quadrant  in  which  the  channel  was  located.  This 
scheme  is  shown  in  Fig  7,2.  In  this  figure  the  wing-to-wing 
(ww)  and  nose-to-tail  will  indicate  the  perpendicular  direc¬ 
tion  of  the  coil  axis.  The  sign  of  the  fields  is  also  shown 
because  these  are  not  according  to  the  usual  coordinate  axis. 

As  mentioned  earlier  in  this  chapter,  the  pitch  and 
banking  angles  of  the  aircraft  can  cause  error  in  direction 
finding.  However,  for  most  of  the  data  collected  during  this 
experiment  the  aircraft  was  at  steady  flight  where  these  types 
of  errors  are  negligible. 

Larger  errors  in  direction  finding  are  caused  due  to 
the  reradiation  of  the  incident  magnetic  field  on  the  air¬ 
craft  structure.  Two  sets  of  crossed  magnetic  loops  installed 
on  the  top  and  bottom  of  the  aircraft  will  yield  different 
bearings.  This  is  because  the  geometry  of  the  aircraft  differs 
in  the  upper  and  lower  portions.  A  detailed  description  of 
this  phenomena  is  given  by  Parker  et  al.,  (1982), 

The  third  error  in  bearing  detection  is  due  to  the 
nonlinear  polarization  of  the  incident  field.  Nonlinearity  in 
polarization  is  caused  due  to  nonhomogeneous  atmospheric  con¬ 
ditions  and  channel  tortuosity.  Bearing  error  due  to  tor¬ 
tuous  channels  near  the  aircraft  can  be  large  and  no  data  are 
available  on  the  amount  of  tortuosity  in  the  channel,  Parker 
et  al.  (1982)  indicates  that  direction  error  in  this  case  is 
given  by 
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Bearing  error  =  tan 
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(7.5) 
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where  L  L  and  L  are  the  components  of  tye  lightning  di- 

A  j  ^ 

pole  in  the  x,  y,  and  z  direction  and  G  is  the  geometrical 
factor 


G  =  (7.6) 

whre  z'  and  z  are  the  altitudes  of  the  dipole  and  aircraft, 
respectively.  D  indicates  the  horizontal  distance  of  the 
dipole  from  the  sensors.  Q  is  a  factor  dependent  on  geo¬ 
metry  and  the  moment  of  the  dipole. 

The  sky  wave  phonomena  as  discussed  by  Gradner  (1980) 
also  introduces  bearing  error  but  in  our  experiment  this 
error  is  also  of  a  negligible  quantity. 

The  conclusion  of  this  discussion  is  that  even  though 
crossed  magnetic  loops  can  be  used  to  det^t  channel  bearing, 
the  results  are  not  accurate.  A  rough  estimate  of  the  chan¬ 
nel  bearing,  however,  can  be  determined  by  this  method. 

Range  Measurement 

In  this  work  tye  electromagnetic  fields  predicted 
by  a  model  are  compared  with  experimental  data  at  corres¬ 
ponding  ranges.  As  calculated  in  Chapter  V,  (equations  5.19 
and  5.3),  the  fields  are  highly  dependent  on  the  range  from 
the  observer  to  the  lightning  channel.  Particularly  in 
the  near  ranges  EM  fields  are  proportional  to  ^  and  ^ 


A  small  error  in  range  detection,  therefore,  compounds  the 
error  in  fields  due  to  exponents  in  R.  Consider,  for  ex¬ 
ample,  the  error  in  range  R  is  5R.  Then  for  the  value  of 
R™  we  get 

(R+SR)"*  =  R®  +  mR®~^  &R 

to  first  approximation  and  the  error  in  R®  is  mR*"  5R. 

Thus  when  m> 1  a  small  error  fiR  in  R  is  increased  in  R®. 

In  actual  data  collection  three  methods  were  used  to 
estimate  the  channel  location  with  reference  to  the  air¬ 
craft.  The  first  method  consists  of  finding  the  bearing  of 
the  channel  as  discussed  above  and  then  correlating  this 
result  with  the  isolated  cloud  source  obtained  from  the 
aircraft  radar,  A  typical  radar  picture  of  the  clouds  is 
shown  in  Fig  7.4.  This  method  of  estimating  the  range  is 
not  accurate  and  errors  of  100%  are  possible  for  large  cloud 
cells.  The  minimum  diameter  of  the  clouds’  in  the  airborne 
data  is  about  4  km. 

The  range  calculated  by  this  method  was  further  com¬ 
pared  with  a  stormscope  display.  The  stormscope  provides 
dots  in  a  cathode-ray-tube  screen  display  whenever  a  large 
discontinuity  occurs  in  the  field.  From  the  crossed  loops 
it  calculates  the  azimuth  to  the  flash  and  the  range  is 
determined  by  using  the  fact  that  the  magnetic  field  is 
inversely  proportional  with  distance.  Under  this  assump- 
tion,  H  =  ^  where  k  is  proportional  to  time  derivative  of 
channel  current  and  considered  to  be  constant.  However,  as 
has  been  discussed  in  the  channel  models,  the  currents  are 


L4 


f 


I.  -r^ 


highly  variable.  Thus,  the  range  predicted  by  the  stormscope 
technique  may  not  be  very  accurate  (Parker  et  al.,  1982). 

From  this  description  it  is  clear  that  the  supporting  in¬ 
formation  for  range  detection  is  not  entirely  accurate  and 
that  its  correlation  with  the  range  from  the  first  method 
does  not  reduce  the  error. 

The  third  method,  which  could  have  provided  a  realistic 
range  of  the  channel,  consisted  of  the  ground  VHF  network  but 
this  data  have  not  been  processed  at  this  time.  Thunder  data 
was  also  recorded  at  the  ground  station  by  using  a  thunder 
microphone.  However  for  most  of  the  data,  the  lightning  flash 
was  located  more  than  20  km  from  the  central  ground  station 
and  thunder  could  not  be  heard. 

Since  all  the  three  methods  discussed  here  failed  to 
provide  a  correct  range  for  the  channel,  we  considered  the 
possibility  of  using  J.A.  Tiller’s  (1975)  statistical  data 
to  estimate  the  ranges.  Tiller  has  suggesj:ed  that  channel 
range  can  be  estimated  from  the  electric  field  measurements 
at  the  ground  level.  For  close  ranges  his  data  is  shown  in 
Fig  7.3.  However,  when  we  checked  his  data  it  was  found  that 
the  statistical  averages  were  highly  erratic.  Table  7.1 
compares  the  suggested  fields  and  ranges  calculated  by  Tiller 
(1975)  with  the  author's  estimates.  Since  these  results  were 
not  consistent,  this  technique  of  estimating  the  range  was 
abandoned.  The  ranges  given  for  the  experimental  data  repre¬ 
sent  the  best  engineering  approximation  possible  based  on 
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Fig  7.3  Range  Approximation  from  Peak  Electric  Field 
.  Value  (J.A.  Tiller,  1975). 


the  airborne  radar  and  the  stormscope  data 


Channel  Height  Problem 

As  is  evident  from  equations  (5. IS)  and  (5.37)/  the  - 
field  at  any  point  depends  upon  the  channel  height  H,  Since 
no  information  was  recorded  about  the  channel  height  it  is 
assumed  that  the  channel  height  was  about  5  km.  This  is  the 
number  used  in  the  euqations  to  compare  the  measured  data  and 
the  predicted  data. 
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VIII.  Comparison  of  Models 

The  linear  and  nonlinear  mathematical  lightning 
return  stroke  models  presented  in  this  thesis  differ  from 
each  other  in  many  aspects.  Also,  the  predicted  results 
from  these  models  are  nonidentical.  To  analyse  models, 
we  first  consider  the  differences  in  their  modeling  pro¬ 
cedure  and  then  we  will  compare  the  results  obtained  from 
them.  It  is  obvious  from  the  description  of  the  models 
that  each  model  has  a  distinct  level  of  sophistication 
which  emerges  from  the  approach  taken  to  describe  the 
lightning  channel  and  the  predicted  EM  fields.  However, 
the  models  presented  in  this  thesis  are  divided  into  three 
categories  on  the  basis  of  the  methods  to  describe  the 
channel  current. 

The  most  simple  and  adequate  approach  in  the 
modeling  of  the  lightning  process  is  to  assume  a  tem¬ 
poral  and  spatial  form  of  the  channel  current.  The 
current  wave  form  and  its  behavior  in  the  assumed  geometry 
of  the  channel  is  derived  from  the  statistical  data  ob¬ 
tained  from  the  natural  lightning.  Since  it  is  difficult 
to  measure  the  current  parameters  in  space,  therefore, 
most  of  the  data  is  collected  at  the  ground  level  and 


from  remote  electromagnetic  field  measurements.  The 
accuracy  of  the  collected  data  depends  on  the  technique 
and  type  of  the  equipment  in  use.  Models  presented  by 
Bruce  &  Golde  (1941),  Dennis  and  Pierce  (1964)  and  Lin 
et  al.  (1980)  were  evolved  by  following  this  technique. 

Lin  et  al.  model  (1980)  has  recently  been  improved  by  Master 
et  al.  (1981)  and  the  newest  model  is  assumed  to  provide 
more  realistic  fields  above  the  ground  level.  The  model 
proposed  by  Gardner  (1980)  is  also  of  this  type  but  differs 
in  the  approach  to  calculate  EM  fields.  The  channel  geom¬ 
etry  in  the  Gardner  (1980)  model  is  of  generalized  form 
and  thus  allows  flexibility  in  calculating  fields  due 
to  tortuous  channels. 

The  second  approach  in  the  modeling  of  lightning 
channel  is  to  consider  the  return  stroke  channel  as  an 
ordinary  discharging  RLC  transmission  line.  The  lumped 
elements,  RLC,  forming  multilooped  transmission  line  as 
shown  in  Fig  4.2  vary  with  height  and  time  and  as  a 
result  of  this,  the  channel  current  becomes  a  function 
of  height  and  time.  The  accuracy  of  this  type  of  model 
depends  on  the  number  of  segments  describing  the  trans¬ 
mission  line.  Uman  (1969)  and  Little  (1978)  have  used 
this  approach  to  model  return  strokes.  However,  Little 
(1978)  has  used  constant  RLC  parameters  in  an  unbranched 
lightning  channel  with  the  assumption  that  current  pulse 
shape  is  identical  at  all  points  along  the  channel.  Due 
to  the  last  assumption,  no  corona  charge  can  be  removed 


from  the  return  stroke  channel  during  the  return  stroke 
propagation  time,  charge  being  only  transferred  from  the 
top  to  the  bottom  of  the  channel  (Lin  et  al.,  1980). 
Parameters  of  lightning  current  pulses  deduced  from  this 
type  of  model  are  in  rough  agreement  with  the  ground  ob¬ 
servation. 

The  third  and  the  most  sophisticated  approach 
in  return  stroke  modeling  follows  the  theory  of  the  de¬ 
velopment  of  an  arc  channel.  Braginskii  (1958)  has 
described  ^he  physic  of  an  arc  channel.  The  arc  channel 
model  is  explained  on  the  basis  of  equations  of  conserva¬ 
tion  of  mass,  momentum  and  energy.  Since  the  channel 
current  and  channel  arc  parameters  interact  with  each 
other  at  all  points  along  the  channel,  therefore,  a  de¬ 
tailed  knowledge  of  parameters,  such  as  ionization  and 
recombination  coefficients,  electrical  and  thermal  con- 
ductives  in  the  channel,  is  necessary  to  ^escribe  the 
model.  Strawe  et  al .  (1980)  has  developed  this  type  of 
lightning  return  stroke.  The  channel  current  is  a  func¬ 
tion  of  environmental  conditions  and  hence,  it  varies 
with  time  and  height. 

Comparison  of  Current  Waveform 

Channel  current  waveform  in  any  return  stroke 


model  is  of  vital  importance  because  the  EM  fields  pi;_- 
dicted  by  the  models  are  strongly  dependent  on  it.  The 


most  important  parameters  in  a  current  waveform  are  the 
risetime,  falltime  and  the  peak  value  of  the  current. 

As  the  physical  conditions  in  a  first  return  stroke  and 
a  subsequent  return  stroke  are  markedly  different,  there¬ 
fore,  the  parameters  which  characterize  the  current  wave¬ 
form  also  differ  in  both  kinds  of  return  strokes.  Fig  2.5 
illustrates  such  differences  as  observed  by  different 
researchers.  One  of  the  earliest  and  most  widely  used 
empirical  approximations  to  the  return  stroke  current 
waveform  at  ground  is  due  to  Bruce  and  Golde  (1941).  The 
waveform  is  graphically  shown  in  Fig  3.1,  whereas,  equa¬ 
tion  3.3  represents  it  in  mathematical  form.  The  round 
peaked  current  waveform  has  average  time  to  crest  and 
to  half  value  of  about  5.6  usee  and  24  microseconds, 
respectively.  The  peak  value  of  the  current  is  21  kA. 

In  Bruce  &  Golde  (1941)  model,  the  current  parameters  in 
first  and  subsequent  strokes  are  considexed  to  be  identi¬ 
cal  and  this  assumption  is  wrong  because  subsequent  stroke 
current  has  a  short  risetime  as  illustrated  in  Fig  2.5. 

Due  to  equation  6.2,  Lin  (1978)  has  found  that  the  return 
stroke  velocity  in  the  Bruce  &  Golde  (1941)  model  is  dis- 
stant  dependent  and  thus,  the  peak  current  in  Bruce  &  Golde 
model  (1941)  increases  sharply  as  the  distance  increases. 

Another  disadvantage  in  Bruce  &  Golde  model  is 
the  infinite  speed  of  information  transformation.  Apart 
from  correcting  this  drawback  in  the  Bruce  &  Golde  model 


(1941),  Dennis  &  Pierce  (1964)  slightly  modified  the 
current  waveform  suggested  in  this  model.  Instead  of 
having  double  exponent  waveform,  Dennis  &  Pierce  (1964) 
suggested  a  triple  exponent  current  waveform  as  in  equa¬ 
tion  2.3.  The  risetime,  falltime  to  50%  of  peak  value 
and  current  peak  value  in  the  first  stroke  due  to  this 
waveform  are  12.8  usee,  52.9  usee  and  21  kA,  respectively. 
But  for  the  subsequent  stroke,  these  values  are  1.0 
usee,  50.1  usee  and  9.8  kA,  respectively.  The  peak 
current  values  agree  with  the  experimental  data,  whereas, 
the  risetime  in  the  waveform  is  excessive.  It  may  be  noted 
that  the  waveform  suggested  either  by  Bruce  and  Golde 
(1941)  or  Dennis  &  Pierce  (1964)  have  poor  representation 
of  the  initial  rise-to-leak  current  when  compared  with 
the  experimental  data. 

The  current  waveforms  derived  from  the  transmission 
line  model  of  Little  (1978)  are  shown  in^fig  4.6a,b,&c.  The 
current  waveform  near  ground  is  sharply  peaked  and  has  a 
very  short  half  value  time.  At  the  ground  level,  the  cur¬ 
rent  risetime  varies  from  1.3  usee  to  3,5  usee  without 
any  significant  frontal  risetime  and  thus,  exceeds  the 
values  found  by  Wiedman  and  Krider:  (1978).  Due  to  the 
recent  studies  of  Weidman  and  Krider  (1978,1980),  the  sharp 
risetime  in  first  and  subsequent  return  strokes  is  iden¬ 
tical  and  of  the  order  of  90  ns.  Also,  the  peak  value 
of  the  current  predicted  by  this  model  exceeds  50  kA  and 
sometimes  it  is  more  than  90  kA;  whereas,  the  measured 


value  of  the  peak  current  is  of  the  order  of  10  kA  (Lin 
et  al.  (1980).  In  view  of  these  current  parameters,  the 
predicted  current  waveform  of  the  Little  model  (1978)  is 
much  different  from  the  measured  waveform. 

Due  to  the  discrepancies  found  in  the  current 
waveform  suggested  in  the  above  mentioned  models,  Lin  et 
al.  (1980)  suggested  a  new  current  waveform,  which  is  com- 
patable  with  the  experimental  data.  As  mentioned  in 
Chapter  III,  the  Line  et  al.  model  (1980)  has  suggested 
three  components  in  the  current  waveform.  In  the  Lin 
et  al.  (1980)  improved  model  (Master  et  al .  1981),  the 
breakdown  pulse  current  in  the  first  stroke  increases  from 
0  to  15  kA  in  5  psec  and  then  it'rises  abruptly  to  30  kA 
in  J.  psec.  The  falltime  to  50%  value  is  approximately 
at  6.9  psec.  Thus,  the  sharp  risetime  of  breakdown  pulse 
is  100  ns.  The  corona  current  waveform  is  similar  to 
Bruce  &  Golde  (1941)  waveform  but,  of  coppse,  with  the 
changed  value  of  a  and  S.  Since  breakdown  pulse  is 
responsible  for  the  peak  value  and  risetime,  thus,  it 
has  removed  the  d..screpancies  in  the  current  forms 
suggested  in  the  earlier  models.  In  Master  et  al.  (1981) 
model,  the  value  of  peak  current  for  the  subsequent  stroke 
is  17.9  kA.  This  includes  the  uniform  current  value,  which 
is  3  kA.  In  the  subseqeunt  stroke,  the  breakdown  pulse 
current  increases  from  0  to  3  kA  in  1.0  psec  and  then  to 
a  peak  value  of  14.9  kA  at  1.1  psec.  It  is  apparent  that 


the  fast  risetimes  in  the- first  and  subsequent  return 
strokes  are  identical.  This  is  exactly  in  accordance 
with  the  findings  of  Wiedman  and  Krider  (1978,1980). 
However,  the  discrepancy  in  Lin  et  al,  (1980)  current 
waveform  is  an  impulse  like  initial  current  peak.  Such 
current  peaks  are  not  observed  in  actual  waveforms.  The 
direct  measurements  of  current  as  shown  in  Fig  2.5  indi¬ 
cate  a  smooth  curved  wave  shape  near  the  peak  values. 

The  curent  waveform  predicted  by  Strawe  et  al. 
model  (1980)  is  shown  in  Fig  8.2.  In  case  of  first 
return  stroke  current  waveform,  there  is  a  slow  rising 
front  of  about  5  psec  duration  and  then  the  current  rea¬ 
ches  to  its  peak  value  of  24  kA.in  a  very  short  time 
(approximately  0.4  psec).  However,  in  the  case  of  a  sub¬ 
sequent  stroke,  there  is  no  slow  rising  front  and  the 
current  reaches  to  its  peak  value  of  22.8  kA  in  about  1.6 
psec.  The  falltime  to  50%  of  peak  value  is,  however, 
very  large  in  this  case  and  this  is  contjrary  to  the  value 
found  from  experimental  data.  Except  for  the  falltime  in 
these  predicted  current  waveforms,  the  parameters  of 
risetime  and  current  peak  values  are  in  good  agreement 
with  the  experimental  data. 

The  nonlinear  model  presented  by  Gardner  (1980) 
does  not  specifically  depend  on  a  current  waveform.  How¬ 
ever  in  the  computation  of  the  EM  field,  the  current  wave 
form  of  Dennis  &  Pierce  (1964)  given  by  equation  2.3  is 
used.  Since  this  waveform  has  a  slow  risetime  and  does 


Fig  8. 


2  Current  Waveform  in  a  First  Stroke  with  Straight 
Verticle  Channel,  Waveforms  are  at  0,  400, 

1200,  1800  and  2200  Meters  (Strawe  et  al.,  1980) 
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not  account  for  the  fast  transition,  therefore,  the  re¬ 
sults  dervied  by  Gardner  (1980)  on  the  basis  of  this  in¬ 
appropriate  current  waveform  are  inaccurate.  The  flex¬ 
ibility  in  the  selection  of  current  waveform  in  Gardner's 
model  promises  for  the  better  results.  Gardner's  model, 
in  fact,  is  more  concerned  about  the  finite  conductivity 
of  the  earth's  surface  and  the  reflection  of  the  radiation 
field  from  the  ionosphere  than  the  current  waveform. 

Table  8-1  shows  a  comparison  of  currents 
waveform  parameters  for  the  models  presented  in  this 
thesis. 

The  behavior  of  the  current  pulse  during  its  up¬ 
wards  propagation  in  the  channel  is  also  a  point  of  inter¬ 
im  est  in  the  return  stroke  modeling.  Except  for  the  models 

presented  by  Lin  et  al.  (1980)  and  Master  et  al  (1981), 
the  velocity  of  porpagation  of  the  current  pulse  de¬ 
creases  with  height.  The  latest  studies  §how  that  for 
the  first-stroke  the  mean  velocity  within  about  1  km  near 

7 

ground  is  9.6x10  m/s,  and  the  subsequent  mean  stroke 

Q 

velocity  is  1.2x10  m/s  (Uman  et  al.,  1982).  Since  the 
mirror  image  effect  as  mentioned  in  Chapter  VII  is  rarely 
observed  in  the  first  stroke  and  almost  never  in  the 
subsequent  stroke, therefore,  it  is  believed  that  current 
peak  amplitudes  decreases  with  height.  The  current  wave¬ 
forms  predicted  by  Strawe  et  al.  (1980)  model  and  shown  in 
Fig  8.1  and  8.2  clearly  substantiates  this  effect.  In  this 
model,  the  peak  value  of  the  current  at  ground  due  to  a 
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first  and  subsequent  stroke  reduce  to  almost  33%  value 
at  the  height  of  1800  meters.  However,  the  attenuation 
factor  for  current  waveform  in  the  Master  et  al.  model 
(1981)  is  a  function  of  height  and  its  value  differs  for 
the  first  and  the  subsequent  stroke.  In  Bruce  and  Golde 
model  (1941),  the  current  is  uniform  at  any  instant  and 
thus  the  peak  value  of  the  current  at  any  altitude  is  the 
same  as  is  obvious  from  Fig  3.1.  This  discrepancy  in  the 
current  waveform  also  exists  in  the  model  of  Dennis  & 
Pierce  (1964).  From  Fig  8,1  it  can  be  observed  that 
peak  value  of  current  in  Little  model  (1978)  attenuates 
with  the  increasing  channel  height.  In  three  kilometers 
trip  of  the  return  stroke,  the  peak  value  has  reduced  to 
30%  of  the  initial  peak  at  ground  level.  This  comparison 
shows  that  Bruce  Si  Golde  model  (1941)  and  Dennis  &  Pierce 
model  (1964)  are  of  a  primitive  nature  and  lack  in  the 
basic  definition  of  the  current  waveform,'  The  current 
waveforms  predicted  by  nonlinear  models  of  Little  (1978) 
and  Strawe  et  al.  (1980)  are  almost  similar  in  shape  but 
differ  in  magnitudes.  The  fall time  to  50%  peak  value  is 
better  in  Little's  model  (1978)  than  that  of  the  Strawe's 
model  (1980).  However,  the  uniform  current  postulated 
in  Lin  et  al,  model  (1980)  is  assumed  to  be  flowing  due 
to  the  leader  current  or  it  may  start  soon  after  the  re¬ 
turn  stroke  initiation.  There  is  no  evidence  given  in 
the  Lin  et  al.  model  that  shows  it  is  different  from  the 
corona  current. 
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Since  Gardner  (1980)  has  used  the  same  current 
waveform  as  suggested  by  Dennis  &  Pierce  (1964),  there¬ 
fore,  it  is  needless  to  say  that  his  model  has  made  any 


significant  contribution  to  improve  the  return  stroke 
model . 

Comparison  of  Predicted  EM  Fields 

As  the  computer  solution  to  the  models,  except 
for  Master  et  al  model  (1981),  were  not  available 
during  this  study,  therefore,  such  models  are  compared 
theoretically  for  their  predicted  EM  fields.  The  par¬ 
ticular  characteristics  of  these  models  has,  however, 
already  been  discussed  during  the  presentation  of  these 
models.  The  primary  object  of  this  section  is,  there¬ 
fore,  to  determine  the  degree  of  correctness  in  their 
predicted  fields. 

Although  Bruce  &  Golde  model  (1941)  was  modi¬ 
fied  by  Dennis  and  Pierce  (1964)  to  avoid  the  require¬ 
ment  of  instantaneous  transfer  of  charge,  the  fields 
predicted  by  these  models  are  alike.  This  is  because 
the  basic  structure  of  both  the  models  is  the  same. 

Lin  (1978)  has  tested  Bruce  &  Golde  model  (1941)  and 
according  to  his  findings,  the  predicted  fields  at 
ground  level  for  the  distances  less  than  5  km  matched 
the  experimental  data.  However,  for  distances  greater 
that  5  km  there  is  a  considerable  mismatch  in  the  pre¬ 
dicted  and  experimental  data.  Since  the  Bruce  8c  Golde 


model  (1941)  and,  for  that  matter,  the  Dennis  &  Pierce.^ 
model  (1964)  has  slow  risetime  in  the  current  pulse, 
radiation  term  does  not  predominate  in  the  first 
20  psec.  Thus,  inductive  and  electrostatic  field 
components  become  stronger ,  Due  to  these  reasons,  the 
Bruce  &  Golde  (1941)  and  Dennis  &  Pierce  (1964)  models 
perform  well  for  short  ranges  but  fail  to  perform  well 
for  distant  observation  points  where  radiated  field  is 
dominant.  From  the  data  of  Table  8.1,  it  is  evident 
that  risetime  for  first  return  stroke  in  Dennis  &  Pierce 
model  (1964)  is  larger  than  that  of  Bruce  and  Golde  .. 
model  (1964).  Thus,  for  unbranched  first  return  stroke, 
the  performance  of  Bruce  &  Golde  model  (1964)  is  com¬ 
paratively  better  than  that  of  Dennis  &  Pierce  model 
(1964). 

Little's  model  (1978)  has  comparatively  sharp 
risetime  in  the  current  pulse.  Therefore,  on  the  basis 
of  above  arguments,  its  performance  is..better  for  long 
distances  and  predicted  data  will  not  match  the  experi¬ 
mental  data  for  short  distances.  Also,  as  seen  in 
Fig  8.1,  the  decay  time  of  the  current  pulse  increases 
with  the  height  and  thus,  magnetic  field  does  not  decay 
for  a  considerable  time.  This  is  probably  true  because 
ten  samples  of  experimental  data  shown  in  Fig  8.3  to 
8.12  indicate  considerable  amount  of  magnetic  field 
up  to  200  psec. 
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In  the  model  presented  by  Gardner  (1980),  the 
current  waveform  is  due  to  Dennis  &  Pierce  (1964)  and 
hence,  even  if  the  earth  is  considered  as  perfect  con¬ 
ductor  and  no  ionosphere  reflections  occur,  the  fields 
predicted  are  those  for  the  Dennis  &  Pierce  model.  Since 
the  ionosphere  reflection  does  not  affect  the  data  for 
lower  ranges,  therefore,  Gardner's  model  is  expected 
to  work  satisfactorily  for  low  ranges.  The  effect  of 
ionospheric  reflection  is  dominant  for  frequencies 
more  than  10  KHz.  High  frequency  components  in 
radiated  field  gets  attenuated  for  distances  greater 
than  the  channel  height  and  thus,  overall  field  value 
drops.  This  conclusion  is  in  harmony  with  the  experi¬ 
mental  data. 

In  Strawe's  model  the  current  pulse  has  a  very 
short  risetime  as  compared  to  other  models  and  the 
radiated  field  is  much  more  dominant  as'compared  to 
any  other  model  but  the  results  produced  by  this  model 
are  not  realistic  (Uman  et  al.,  1982). 

The  Masters'  model  (1981)  computer  orogram, 
which  was  made  available  to  the  author,  was  used  to 
compute  the  electric  and  magnetic  fields  at  the  corresponding 
altitude  and  distances  where  experimental  data  were 
collected  by  the  Air  Force  Flight  Dynamics  Laboratory. 

This  data  are  presented  in  Figs  8.3  to  8.10.  In  the 
empirical  data  the  vertical  component  of  electric 
field  (E^)  was  recorded  through  two  sensors.  These 


sensors  were  located  at  aft  upper  fuselage  (AUF)  and 
aft  lower  fuselage  (ALF)  of  the  aircraft.  For  the  near 
ranges,  that  is,  from  27  km  to  11  km. distances,  the  data 
is  only  available  due  to  ALF  sensor.  The  author  was  told 
that  there  is  a  ratio  of  2.6  in  ALF  and  AUF  sensors  data 
but  in  the  available  data  no  such  relation  exists.  Due 
to  this  reason,  only  ALF  sensor  data  are  used  for  the 
analysis.  In  this  data,  electric  field  units  are  in 
volts/meter.  Magnetic  field  data  were  recorded  due  to 
two  orthogonal  sensors  as  discussed  in  Chapter  VII, 
However,  the  peak  magnetic  field  is  the  square  root  of 
the  sum  of  the  squares  of  fields  sensed  by  individual 
sensors.  Risetime  for  magnetic  field  is  calculated 
from  the  best  available  plot.  The  magnetic  field  is 
plotted  in  ampere/meters. 

The  predicted  fields  due  to  Master  et  al,  model 
(1981)  are  shown  in  Fig  8.13  to  8.15,  In  these  figures, 
the  left  side  column  is  for  the  ranges  2.7  km,  7  km  and 
11  km  where  lowest  range  is  on  the  top.  The  right 
side  column  corresponds  to  4  km,  10.2  km  and  30  km 
ranges  with  the  lowest  range  on  the  top.  The  bottom 
plot  is  for  35  km  range.  This  scheme  of  data  presenta¬ 
tion  is  common  in  Figs  8,13  to  8,21,  The  Figs  8,16  to 
8,18  are  computer  plots  for  the  Master  et  al.  model 
with  velocity  variation  as  suggested  in  equation  6.3. 
Figus  8.19  to  8.21  are  the  results  when  current  wave¬ 
form  in  the  Masters’  model  is  modified.  In  predicted 
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data  the  electric  field  is  in  volt /meter,  whereas,  the 
magnetic  field  is  shown  in  amperes/meter , 

The  comparison  of  empirical  data  and  predicted 
data  is  based  on  three  factors,  which  are,  risetime, 
initial  peak  value  of  fields  and  the  time  at  which  the 
electric  and  magnetic  fields  cross  the  base  line.  Most 
of  this  study  is  restricted  to  first  return  stroke  data 
because  the  empirical  data  for  subsequent  return  strokes 
were  not  available. 

The  risetime  in  the  data  is  considered  for  the  por¬ 
tion  where  fast  transition  occurs.  For  the  Master  et 
a 1.  model  (1981),  this  transition  time  is  fixed  and  equal 
to  0.1  usee.  In  the  empirical  data  as  shown  in  Table 
8.2,  the  fast  transition  period  is  slightly  more  than 
1  usee  for  the  ranges  less  than  30  km.  However,  for 
the  ranges  of  30  and  35  km,  the  fast  transition  period 
is  four  times  that  predicted  by  Masteu.^t  al.  model. 

As  the  author's  of  Master  et  al.  model  (1981)  claim 
that  the  model  works  up  to  10  km  ranges,  the  risetime 
in  predicted  data  and  in  the  empirical  data  are  in 
excellent  agreement.  The  excess  value  of  risetime  in 
empirical  data  can  be  attributed  to  measurement  error. 

As  the  risetime  in  Master  et  al,  model  (1981)  was  not 
modified,  therefore,  the  above  result  holds  good  for 
the  modified  models. 


134 


In  empirical  results  electric  field  (E^)  never 
crosses  the  base  line.  Rather,  it  stays  much  above  the 
initial  value.  However,  the  electric  field  predicted 
by  Master  et  al,  model  crosses  the  base  line  for  the 
ranges  less  than  7  km  and  greater  than  11  km.  Thus, 
the  performance  of  thisr  model  is  restricted  to  limited 
ranges  in  respect  of  base  line  crossing.  In  fact, 
as  can  be  seen  from  equation  5.36,  electric  field  can 
cross  the  base  line  when  ^  changes  its  sign  from  +ve 
to  -ve  and  electrostatic  and  induced  fields  are  less 
dominant.  Therefore,  it  seems  that  radiation  field 
becomes  more  dominant  in  the  low  ranges  as  compared  to 
the  experimental  data.  Particularly  this  phenomena 
cocurs  due  to  the  fast  decaying  of  the  corona  pulse. 

In  Fig  8.13,  the  field  due  to  corona  current  has 
crossed  the  base  line  for  the  range  of  11  km,  whereas, 
in  Fig  8.16  it  is  seen  that  due  to  veltJcity  variation, 
this  tendancy  of  crossing  the  base  line  is  reduced  and 
the  peak  value  of  the  field  has  increased.  This  fact 
of  crossing  the  base  line  due  to  corona  current  was 
further  amplified  when  the  corona  current  shape  in 
the  Master  et  al.  model  (1981)  was  modified.  The  results 
of  modified  current  waveform  are  shown  in  Fig  8.19  to 
8.21,  With  the  fast  decaying  corona  current  as  shown 
in  Fig  8.22.  Although  with  this  modification  the 
tendancy  of  crossing  the  base  line  has  increased  but 
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the  peak  value  of  the  electric  field  has  increased  and 
it  is  better  than  the  value  of  electric  field  predicted 
by  Master  et  al.  model  (1981).  The  data  due  to  modified 
velocity  and  current  waveform  are  shown  in  Tableds  8.4 
and  8.5.  It  can  be  inferred  from  these  results  that 
simultaneous  implementation  of  these  modifications  will 
improve  the  peak  field  value  and  base  line  cross 
over  problem.  A  comparison  of  data  in  Tables  8.2  and 
8.3  indicates  that  the  peak  value  of  predicted  fields 
by  the  Master  et  al.  (1981)  model  is  much  less  than  the 
experimental  data.  Also,  the  value  of  I^|  is  almost 
ten  times  less  than  that  for  the  empirical  data  in  the 
near  ranges.  However,  as  the  range  increases,  the 

it 

difference  in  the  peak  values  of  predicted  and  experi¬ 
mental  decreases.  Also,  the  time  derivative  of  the 
electric  field  becomes  almost  identical  to  that  of  the 
empirical  data. 

Another  difference  in  the  data  of  Table  8.2  and 
8.3  is  that  peak  value  of  in  Table  8.3  increases 

from  2.7  km  range  to  7  km  range,  and  then  it  starts  de¬ 
creasing,  This  behavior  of  is  contrary  to 

the  results  from  Table  8,2.  It  is  not  understood 
clearly  why  this  phenomena  occurs  in  Master  et  al. 
model  (1981).  Most  probably  this  is  because  of  the 
image  current  and  channel  height.  Since  after  5  km  the 

"T-  distances  from  the  channel  image  start  becoming  com¬ 

parable  with  the  direct  distances  from  the  channel  and, 
therefore,  E-field  starts  reducing  with  the  distance. 
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TABLE  8.1 

COMPARISON  OF  CURRENT  WAVEFORM  PARAMETERS 
IN  VARIOUS  MODELS 


Model  Presented  First  Return  Stroke 

"t - n — I - 

rip 
lisec  usee  kA 


Subsequent  Return  Stroke 

tl  tZ  x~ 

r  I  p 

usee  usee  kA 


Bruee  &  Golde 
(1941) 


5.6  24  21 


Dennis  &  Pieree  12.8  52.9  21 

(1964) 


1.0  50.1  9.8 


Master  et  al. 
(1981) 


5.1  6.9  35 


3.8  17.9 


Little  (1978)  1.3-  40.0  50-90 

3.5 


Gardner  (1980) 


Same  as  for  Dennis  &  Pieree  (1964) 


Strawe  et  al.  1.6  >100  24 

(1980) 


>100  11.8 


t^  =  risetime 


tj  =  fall  time  to  50%  of  peak  value 
i„  =  peak  value  of  the  ehannel  eurrent 
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TABLE  8.3 


DATA  FROM  THE  FIRST  RETURN  STROKE 
(Master  et  al.  Model,  1981) 


D 

km 

E 

km 

re  rb 
ysec 

!e  I 

'  max ' 
V/m 

l®n.axl 

A/m 

dE 

dt 

V/m-sec 

,dB| 

•at' 

A/m-sec 

2.7 

5.1 

0.1 

13.7 

.142 

1.17x10® 

0.69x10^ 

4.0 

5.1 

0.1 

31.74 

.16 

1.8x10® 

0.97x10® 

7.0 

5.1 

0.1 

45.2 

.15 

2.39x10® 

0.68x10® 

10.2 

5.1 

0,1 

42.0 

.12 

2.15x10® 

0.55x10® 

11.0 

5.1 

0.1 

40.12 

.073 

2.1x10® 

0.244x10® 

30.0 

4.0 

0.1 

17.24 

.0447 

0.81x10® 

0.225x10® 

35,0 

4.0 

0.1 

14.3 

.0397 

0.74x10® 

0.202x10® 

Symbols  are  same  as  in  Table  8.2 


TABLE  8-4 


I 


DATA  FOR  FIRST  RETURN  STROKE  WITH 
: MODIFIED  PROPAGATION  VELOCITY 


1  dEi 

1  db  1 

D 

z 

t  4t  , 
re  rb 

1 ^max  ^ 

B 

'  max ' 

latl 

Idtl 

km 

km 

usee 

V/M 

A/M 

V/M-sec 

A/m-sec 

2,7 

5.1 

0.1 

25 

0.216 

2.5x10° 

10.4x10 

4.0 

5.1 

0.1 

44.4 

0.23 

2.5x10® 

12.2x10 

7.0 

5.1 

0.1 

63.5 

0.216 

3.3x10® 

11.8x10 

10.2 

5.1 

0.1 

50.8 

0.176 

2.3x10® 

9.5x10® 

11.0 

5.1 

0.1 

56.3 

0.167 

3.1x10® 

9.07x10 

30.0 

4.0 

0.1 

24.56 

0.064 

1.31x10® 

3.33x10 

35.0 

4.0 

0.1 

21.05 

0.056 

1.14x10® 

3.08x10 

Symbols 

are  the 

same  as 

in  Table 

8r2 

TABLE  8-5 


DATA  FOR  FIRST  RETURN  STROKE  WITH 
MODIFIED  CURRENT  PULSE  WAVEFORM 


D 

km 

Z 

km 

psec 

IE  II 

•  max  >  1 

V/m 

®maxl 

A/m 

1^1 

Idtl 

V/m-sec 

1^1 

Idtl 

A/m-sec 

2.7 

5.1 

0.1 

41.22 

.35 

3.42x10® 

18.15x10' 

4.0 

5.1 

0.1 

69.84 

.38 

3.49x10® 

21.5x10® 

7.0 

5.1 

0.1 

88.88 

.307 

4.76x10® 

16.75x10 

10.2 

5.1 

0.1 

86.792 

.237 

4.9x10® 

13.26x10 

11.0 

5.1 

0.1 

78.07 

.244 

4.2x10® 

12.24x10 

30.0 

4.0 

0.1 

33.04 

.086 

1.7x10® 

3.91x10® 

35.0 

4.0 

0.1 

28.94 

.073 

1.6x10® 

1.03x10® 

Symbols 

have 

the  same  meaning  as  in 

Table 

8-2 
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Fig  8.3  Electric  and  Magnetic  Fields  for  a  First  Stroke 
At  a  Distance  of  2.7  km  and  Altitude  5,1  km. 
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Fig  8.4  Electric  and  Magnetic  Fields  for  the  First  Stroke 
at  a  Distance  of  4  km  and  Altitude  5,1  km 


Fig  8.6  Electric-  and  Magnetic  Fields  for  the  First 
Return  Stroke  at  a  Distance  of  10.2  km  and 
Altitude  5.1  km 
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Fig  8,9  Electric  and  Magnetic  Fields  for  a  First  Return 
Stroke  at  a  Distance  of  30  km  and  Altitude  4  km 
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Fig  8.10  Electric  arid  Magnetic  Fields  for  a  First  Stroke 
at  a  Distance  of  30  km  and  Altitude  4  km. 
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Fig  8,13  Verticle  E-Field  due  to  Master  et  al 
(1981).  Range  (km);  (a)  2.7;  (b)  4.6 
(d)  10.2;  (e)  11.0;  (f)  30.0;  (g)  35. 
(km):  5  for  a-e,  4  for  g-h. 


-  Model 
;  (c)  7.0 
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IX.  Conclusions  and  Recommendations 

Conclusions 

Master  et  al,  model  (1981)  presented  In  this 
thesis  predicts  fast  rlsetlme  correctly.  However,  the 
peak  value  of  EM  field  Is  much  less  than  empirical  value, 

EM  field  In  actual  data  never  crosses  the  base  line  but 
the  fields  predicted  by  Master  et  al.  model  (1981)  very 
often  cross  the  base  line.  The  modified  velocity  for 
current  pulse  propagation  has  significantly  reduced  the 
tendency  of  crossing  the  base  line.  With  this  modified 
velocity,  the  peak  value  of  the  EM  fields  has  also  Improved 

Modified  current  waveform  has  Increased  the  peak 
value  of  EM  fields  very  much  but  with  this  the  tendency 
of  crossing  the  base  line  has  been  advanced  In  time  scale. 
It  Is  hoped  that  the  combined  effect  of  modified  velocity 
and  current  waveform  will  Improve  the  predicted  values  of 
EM  field. 

Since  no  significant  amount  of  empirical  data  for 
subsequent  strokes  was  available,  therefore,  this  aspect 
of  Master  et  al,  model  (1981)  could  not  be  studied.  Ve¬ 
locity  variation  In  subsequent  strokes  Is  very  nominal. 
However,  from  the  comparison  of  different  models  and 
waveforms  measured  by  different  researchers.  It  is  felt 


that  current  waveform  modification  in  Master  et  al,  model 
is  necessary.  Particularly,  the  peak  value  of  EM  fields 
predicted  by  Master's  model  (1981)  are  significantly  low 
in  space. 

Strawe's  model  (1979)  has  very  good  risetime  in 
its  predicted  current  waveforms.  Also,  as  the  current 
waveform  in  this  model  varies  with  time,  therefore,  this 
model  has  great  potential  to  provide  realistic  results. 

At  present,  this  model  needs  to  be  refined.  Gardner 
model  (1980)  presently  does  not  propose  a  better  solu¬ 
tion  for  the  problems  in  existing  models.  The  concept 
presented  in  this  model  can  be  used  only  for  the  finer 
details  in  a  successful  model. 

Recommendat ions 

In  the  detailed  study  of  Master  et  al,  model  more 
data  samples  should  be  collected  and  compared  to  the  pre¬ 
dicted  results.  Since  this  model  is  primarily  meant 
for  subsequent  strokes,  empirical  data  for  subsequent 
strikes  is  required  to  correctly  assess  the  performance 
of  this  model.  For  the  first  return  stroke,  current 
waveform  is  required  to  be  modified  in  order  to  improve 
the  peak  EM  fields  and  also  to  reduce  the  base  line  cross 
over  tendency. 

Strawe  et  al.  model  (1980)  computer  program  should 
be  modified  for  the  channel  diameter.  The  channel  dia¬ 
ls  2 


meter  reduces  with  time  but  it  can  be  considered  as  a 
function  of  height  as  well.  Thus,  this  dual  variation  in 
channel  diameter  with  height  and  time  is  hoped  to  produce 
better  results  or  at  least  it  will  provide  an  opportunity 
to  assess  the  behavior  of  predicted  EM  field  with  channel 
diameter  variation. 

Gardner's  model  may  be  tested  with  current  waveform 
of  Master  et  al.  model.  The  predicted  result  once  compared 
with  experimental  data  will  allow  the  future  researchers 
to  assess  the  credibility  of  this  model. 
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